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Wave  action  density,  where  A  =  E/a 
phase  speed  for  surface  waves 
group  velocity  of  Eq  (6)  for  surface  waves 
fetch 

local  energy  density 
wind  directional  function  of  Eq  (13) 
energy  spectral  density  of  Eq  (12) 
acceleration  due  to  gravity 
wavenumber  vector  of  magnitude  k 
slope,  see  Eq  (98) 

apparent  or  observed  frequency  of  Eq  (2) 
spectral  action  density,  where  N(k)  =  F(k)/a 
source  term  of  Eq  (19),  see  also  Eq  (20)  ■*  ( 
radiation  stress  as  found  in  Eq  (7) 
time 

ratio  of  surface  tension  to  water  density 
friction  velocity 

displacement  spectral  density  of  Eq  (10) 
,nondimensional  fetch 
location  vector 

;  distance  above  the  mean  sea  surface 


i 


a  ■  .  r|  see  Eq  (20) 

P*  ■  •  ^9  (21) 

*  •%  j  #  I  ^ 

y.  ,  Phillips  constant  of  Eq  (14), 

see  Eq  (22) 


(30) 


A(k,k’) 


attenuation  coefficient  in  presence  of  dense  slick 
attenuation  coefficient  for  clean  surface 
see  Eq  (23) 

surface  displaceisent  about  the  mean  surface 
wavelength,  where  k  =  2n/\ 
kinetic  viscosity 
density 

intrinsic  frequency  whose  k  dependence  gives  the  dispersion 
relation 

displacement  frequency  density  of  Eq  (11a) 
displacement  spectral  density  of  Eq  (11b) 


I .  INTRODUCTION 


Two  inportant  naval  system  problems,  acoustic  torpedo  surface 
reverberation  in  shallow  water  and  remote  radar  sea  surface  imaging,  have  a 
common  physical  interface,  ie,  the  sea  surface.  The  proper  and  adequate 
modeling  of  the  sea  surface  determines  our  ability  to  model  and  predict  the 
behavior  of  homing  torpedoes  in  shallow  water  and  the  imaging  of  remote  radar 
sensors.  Therefore,  this  report  examines  the  state  of  the  art  in 
ocean-atmosphere  interface  modeling. 

A  survey  of  recent  oceanographic  literature  showed  that  there  are 
relatively  few  fundamental  mathematical  relations  for  modeling  the  sea 
surface.  These  relations  are  presented  in  Section  II  and  applied  in  Section 
III.  In  addition,  several  useful  theoretical  results  and  einpirical  relations 
are  included  throughout  this  report.  As  shown  throughout  Section  III,  many 
models  begin  with  the  appropriate  fonss  of  the  kinematical  equation  (Section 
II. A.)  and  the  action  conservation  principle  (Section  II.B.).  Sometimes  the 
radiative  stress  equation  (Section  II. C)  is  used  instead  of  the  action 
conservation  principle.  A  much  older  and  less  prolific  relation,  the 
Bernoulli  equation,  is  discussed  in  one  of  its  useful  forms  in  Section  II. D. 

The  waves  found  on  the  surface  of  the  sea  are  almost  always  random  in  the 
sense  that  the  detailed  configuration  of  the  sea  surface  varies  in  an 
irregular  manner  in  both  space  and  time.  Therefore,  only  the  various 
statistical  measures  of  the  surface  motion  can  be  regarded  as  significant 
observationally  or  predictable  theoretically.  Consequently,  Sections  II. E  and 
II. F  discuss  various  spectral  relations  (and  their  directional  behavior)  and 
concepts.  Equation  (18),  the  conservation  of  action  spectral  density,  often 
proves  to  be  a  useful  statistical  relation  for  sea  surface  modeling  purposes. 
However,  the  radiative  transfer  equation,  Eq  (19),  promises  to  be  the  most 


fruitful  statistical  relation  for  modeling  the  sea  surface  and  is  discussed  in 
Section  II. F. 

The  state  of  the  art  in  practical  ocean  surface  modeling  is  treated  in 
Section  III  and  represents  the  main  results  of  this  survey.  The  "mean  JONSWAP 
spectrum"  (Section  III.B)  is  a  well-accepted  relation  and  provides  a  major 
empirical  benchmark  for  many  theoretical  developments.  This  study  was 
designed  to  treat  the  effect  of  four  phenomena  on  the  sea  surface  behavior. 
These  four  phenomena  are  wind  velocity,  surface  current,  swell,  and  internal 
waves.  The  directional  behavior  of  the  wind  is  discussed  in  Section  II. E  and 
the  effects  of  its  magnitude  are  shown  in  Section  III.B.  Four  significant 
properties  of  a  variable  surface  current  are  presented  in  Section  III.C. 

These  are  (1)  upwelling,  (2)  lateral  spreading,  (3)  frequency-dependent 
attenuation,  and  (A)  refraction  of  surface  waves  by  a  shear  current.  Section 
III.D  discusses  the  interaction  of  swell  and  short  gravity  waves  and  gives 
some  special  consideration  to  imaging  by  a  remote  radar  sensor.  The  most 
thoroughly  developed  and  verified  theory  concerns  the  interaction  of  surface 
and  internal  waves.  This  is  discussed  in  Section  III.E.  Finally,  Section 
III.F  considers  wave  attenuation. 

In  Section  IV  it  is  tentatively  concluded  that  the  state  of  the  art  in 
sea  surface  modeling  is  sufficiently  developed  to  permit  the  two  important 
naval  system  problems  specified  to  be  addressed.  Therefore,  Section  IV 
proceeds  to  suggest  a  possible  method  of  developing  a  software  model  for  this 
purpose.  The  appendix  examines  imaging  for  remote  radar  sensing. 


II.  FUNDAMENTAL  MATHEMATICAL  RELATIONS 


A.  KINEMATICAL  EQUATION 

The  kineaatical  equation  (Ref  1)  gives  the  kinematical  conservation  of 
the  density  of  waves  in  the  form 

8k 

=  0  (1) 

as  given  by  Ref  2  and  3,  where  k  is  the  wavenumber  vector  and 

n  =  a  +  k*U  (2) 

is  the  apparent  or  observed  frequency  of  waves  passing  a  fixed  point  when  the 
medium  itself  is  moving  with  velocity  U(x,t),  which  is  a  slowly  varying 
function  of  position  x  and  time  t,  and  a  =  a(k)  is  the  intrinsic  frequency 
whose  k  dependence  represents  the  dispersion  relation.  The  distribution  of 
the  local  wavenumber  k  in  space  obeys 

V  X  k  =  0.  (3) 

In  deep  water 

a2(k)  =  gk  +  Tk®,  (4) 

where  g  is  the  acceleration  due  to  gravity,  and  T  is  the  ratio  of  the  surface 
tension  to  the  water  density  (Ref  1,  p  37-38).  The  first  right-hand  term  in 

1.  O.M.  Phillips,  The  Dynamics  of  the  Upper  Ocean,  Cambridge  University 
Press,  1977. 

2.  F.  Unsell,  "Steady  wave  patterns  on  a  non-uniform  fluid  flow," 

J.  Fluid  Mech. .  9,  1960,  333-346. 

3.  G.B.  Whitham,  "A  note  on  group  velocity,"  J.  Fluid  Mech.,  9,  1960, 
347-352. 
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£q  (4)  represents  gravity  waves,  while  the  second  represents  capillary  waves; 
gravity  waves  dominate  for  wavelengths  X  >  10  cm,  where  k  =  2n/\. 

B.  ACTION  CONSERVATION  PRINCIPLE 

The  action  conservation  principle  (established  by  Bretherton  and  Garrett 
(Ref  4)  and  elaborated  further  by  Bretherton  (Ref  5);  see  also  Ref  1, 
p  26-27]  is 

+  V-I(U  +  c^)«A]  =  0,  (5) 

ot  -  -g 

where  the  wave  action  density  A  =  E/a,  E  is  the  local  energy  density,  and 

c  =  V.a(k)  (6) 

-g  k 

is  the  group  velocity  (Ref  1,  p  25).  Equation  (5)  is  valid  generally  for 
nondissipative,  short,  small -amplitude  waves  in  a  moving  medium.  It  states 
that  the  local  rate  of  change  of  A  is  balanced  by  the  convergence  of  the  flux 
of  A  that  flows  relative  to  the  moving  medium  with  velocity  c^.  It  can  be 
shown  that  the  phase  speed  c  for  a  surface  gravity  wave  obeys  c^  =  %c. 

C.  THE  RADIATION  STRESS  EQUATION 

If  the  fluctuating  motion  consists  entirely  of  a  wave  train  in  which  the 
energy  dissipation  is  negligible,  the  energy  balance  is  given  by 

4.  F.P.  Bretherton  and  C.J.R.  Garrett,  "Wavetrains  in  inhomogeneous  moving 
media,”  Proc.  Roy.  Soc.,  A,  302,  1969,  529. 

5.  F.P.  Bretherton,  "The  general  linearized  theory  of  wave  propagation,"  in 
Mathematical  Problems  in  the  Geophysical  Sciences,  vol.  1,  61-102,  Amer. 

Math.  Soc.,  Providence,  RI,  1971 


0& 


(7) 


a 

+  s^{E[U  +  (c  )  1}  +  S  «  = 

at  ax  *  ^  a  '  a^  ax^ 

a  Of 


0 


(Ref  1,  p  63-68,  and  Ref  6-8),  where  a,P  =  1,2  and  the  radiation  stress 

represents  the  excess  momentuio  flux  associated  with  the  wave  motion  alone. 

Phillips  (Ref  1,  p  68)  shows  that  S  .  s  S  =  ^  in  deep  water. 

“P 


D.  THE  BERNOULLI  EQUATION 

Starting  with  the  momentum  equation  in  the  form 

^  +  U)  X  u  +  V(p/p  +  +  gz)  =  vV^u  (8) 

at 

where  p  is  the  pressure,  p  is  the  density,  V  is  the  kinetic  viscosity,  z  is 
distance  above  the  mean  sea  surface,  w  =  V  x  u  is  the  vorticity,  and 
expressing  the  particle  velocity  u  in  terms  of  the  velocity  potential  or  u  = 
gives  (Ref  1,  p  19-20)  the  Bernoulli  equation  for  irrotational  (or 
potential)  flow  with  w  =  0  and  v  =  0 

^  H  +  %u2  +  gz  =  f(t),  (9) 

where  f(t)  is  an  arbitrary  function  of  time  determined  by  the  pressures 
imposed  at  the  boundaries  of  motion. 


6.  M.S.  Longuet-Higgins  and  R.W.  Stewart,  "Changes  in  the  form  of  short 
gravity  waves  on  long  waves  in  tidal  currents,"  J.  Fluid  Mech. ,  8,  1960, 
565-583. 

7.  M.S.  Longuet-Higgins  and  R.W.  Stewart,  "The  changes  in  amplitude  of  short 
gravity  waves  on  steady  non-uniform  currents,"  J.  Fluid  Mech.,  10,  1961, 
529-549. 


8.  G.B.  Whitham,  "Mass  ammentum  and  energy  flux  in  water  waves,"  J.  Fluid 
Mech. .  12,  1962,  135-147. 


E.  SOME  SPECTRAL  AND  DIRECTIONAL  RELATIONS 


The  mean  square  surface  displacement  is  represented  spectrally  by 

P  =  J  X(k,n)dkdn  (10) 

00 

=  /  (J)(n)dn  =  /  »l<(k)dk  (11) 

o 

so  that 

E  =  /  F(k)dk  (12) 

=  P&U  =  Pg  /  i|<(k)dk. 

From  the  directional  material  in  Barnett  and  Kenyon  (Ref  9,  p  688-689, 
713-714)  and  the  "saturation  range"  material  in  Phillips  (Ref  1,  p  102-104, 
140-149),  the  following  useful  relations  can  be  obtained.  Waves  leaving  a 
storm  area  can  have  a  broad,  say  ±45®,  directional  distributioi  relative  to 
the  mean  wind  direction.  Represent  the  directional  spectrum  (|)i(n,0)  by 

4.i(n,0)  =  <t)(n)f(0)  ,  (13) 

where  0  =  0®  in  the  mean  wind  direction.  For  the  saturation  range  Phillips 
(Ref  1,  p  102-104,  149)  gives 

<|)(n)  =  Yg^n  ^  ,  n^  «  n  «  2g/u^.  ,  (14) 

where  n^  is  the  spectrum  peak  apparent  frequency,  u^.  is  the  friction  velocity, 
and  Y  is  a  "constant"  whose  dimensionless  fetch  dependence  is  given  by  Eq  (30) 
and 


9.  T.P.  Barnett  and  K.E.  Kenyon,  "Recent  advances  in  the  study  of  wind 
waves,"  Rep.  Prog.  Phys.,  38,  1975,  667-729. 


(15) 


i|»(k)  =  Dk'S(0),  k^  «  k  «  2g/u|  , 

where  k  is  the  spectrum  peak  wavenumber  magnitude,  0  specifies  the  direction 
o 

of  k  relative  to  the  mean  wind  velocity,  and  D  ~  (Ref  1,  p  148-150).  For 
an  extension  of  Eq  (14)  to  higher  frequencies  see  Ref  10.  For  a  critique  of 
Eq  (14)  see  Ref  11. 

The  distribution  of  wave  energy  in  k-space  for  a  fully  aroused  sea  was 
first  estimated  with  some  accuracy  during  the  SWOP  project  (Ref  9,  Section 
4.2)  by  Cote,  et  al  (Ref  12).  The  rough  sea  under  the  action  of  a  steady  wind 
gave 

f(0)  =  1  +  aje  ^cos  0  +  32^  ^cos  40  ,  (16) 

4 

where  the  a^  are  constants  and  b  =  %(oU/g)  .  Basically  this  equation  says 
that  the  waves  are  distributed  as  cos^6  near  the  spectral  peak,  with  the 
beamwidth  broadening  at  higher  frequencies.  A  substantial  advance  in 
estimating  f(0)  was  made  by  Longuet-Higgins  et  al  (Ref  13)  and  Cartwright  and 
Smith  (Ref  14): 

f(0)  =  cos*^^^,  -  n  <  0  <  7t  and  s  =  s(n).  (17) 

10.  H.  Mitsuyasu,  "Measurement  of  the  high-frequency  spectrum  of  ocean 
surface  waves,"  J .  Phys .  Ocean . ,  7,  1977,  882-891. 

11.  S.A.  Kitaigorodskii,  "The  statistical  characteristics  of  wind-generated 
short  gravity  waves,"  in  Spaceborne  Synthetic  Aperture  Radar  for  Oceanography, 
R.C.  Beal,  et  al.,  eds.,  The  Johns  Hopkins  University  Press,  1981. 

12.  C.S.  Cote,  et  al..  Meteorological  Papers,  2(6),  W.J.  Pierson,  Jr.,  ed. , 
1960. 

13.  M.S.  Longuet-Higgins,  D.E.  Cartwright,  and  N.D.  Smith,  Ocean  Wave 
Spectra.  Prentice-Hall,  1963,  111-136. 

14.  D.E.  Cartwright  and  N.D.  Smith,  Buoy  Technology,  Marine  Technology 
Society,  Washington,  D.C.,  1964,  112-121. 
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The  values  of  s  ranged  from  1  for  high  frequencies  to  10  at  low  frequencies. 
Ewing  (Ref  15)  confirmed  these  results.  Phillips  (Ref  16)  suggests  f(6)  = 
cos^B  or  cos^B  as  a  rough  rule  of  thumb. 


A  further  useful  spectral  equation  is 

I^N(k)  +  (Cg  +  U)-VN(k)  =  0,  (18) 

where  the  spectral  energy  density  is  related  to  the  spectral  action  density  by 
N(k)  =  F(k)/a  (see  Ref  1,  p  181).  Equation  (18)  expresses  the  conservation  of 
action  spectral  density  N(k,x,t)  of  short  waves. 


F.  THE  RADIATIVE  TRANSFER  EQUATION 

The  radiative  transfer  equation,  Eq  (19),  describes  the  energy  balance  of 
the  wind  wave  field  in  terms  of  the  energy  spectrum  F(k)  of  surface  gravity 
waves.  This  equation,  which  summarizes  all  the  various  physical  processes 
which  can  change  the  wave  energy,  is  given  by 


dt  “  at  8t  ax.  at  ak."  ^i  ’ 

1  1  TTr 


(19) 


where  the  spectrum  F(k,x,t)  is  locally  a  function  of  k  but  is  allowed  to  vary 
slowly  as  a  function  of  x  and  t.  The  characteristic  equations  are 


ax . 
_ 1 

at 


aa 

ak.  ‘'gi 


(6a) 


15.  J.A.  Ewing,  J.  Marine  Res.,  27 ,  1969,  163-171. 

16.  O.M.  Phillips,  "The  structure  of  short  gravity  waves  on  the  ocean 
surface,"  in  Spaceborne  Synthetic  Aperture  Radar  for  Oceanography,  R.C.  Beal, 
et  al,  eds..  The  Johns  Hopkins  University  Press,  1981. 


and 


dk. 
_ 1 

at 


an 

ax. 

1 


(la) 


and  are  equivalent  to  Hamilton's  equation  for  a  particle.  Barnett  and  Kenyon 
(Ref  9,  p  678-684,  698-704,  718-722)  survey  the  state  of  the  art  of  applying 
Eq  (19).  Hughes  (Ref  17)  applies  Eq  (19)  to  the  interaction  of  internal  waves 
and  surface  wind  waves. 

Hasselmann  (Ref  18)  and  Bai-..c:  t  and  Kenyon  (Ref  9,  p  678-682)  give  some 


forms  for  the  source  functions  S^,  nasiely 

Sj  =  a  (20) 

Sa  =  pF(k)  (21) 

53  =  F(k)  /  Y(k,k')F(k')dk’  (22) 

54  =  -6F(k)  +  /  £(k,k')F(k’)dk’  (23) 

85=/  lTiF(k’)F(k-k'-k")  -  T2F(k)F(k’ )F(k")Jdk’dk"  .  (24) 


The  Phillips  linear  growth  function  represents  the  constant  energy  transfer 
to  the  wave  field  through  turbulent  atmospheric  pressure  fluctuations.  The 
Niles  exponential  growth  function  Sa  represents  the  increasing  transfer  of 
energy  to  the  wave  field  due  to  an  instability  in  the  coupling  between  the 
wave  field  and  the  mean  boundary  layer  flow  in  air.  The  term  S3  is  a 
nonlinear  correction  to  83.  The  term  S4  represents  the  energy  transfer  due  to 
the  interaction  between  waves  and  turbulence  in  the  atmosphere.  The  term  85 


17.  B.A.  Hughes,  "The  effect  of  internal  waves  on  surface  wind  waves.  2: 
Theoretical  analysis,"  J.  Geophys.  Res.,  83,  1978,  455-465. 

18.  H.  Hasselmann,  Basic  Developments  in  Fluid  Dynamics.  Vol.  2,  M.  Holt, 
ed..  Academic  Press,  1968,  117-182. 
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represents  the  energy  transfer  among  the  various  wavenumber  components  due  to 
weak  nonlinear  wave-wave  interactions;  85  is  further  treated  in  Ref  9, 
p  682-684,  and  is  discussed  in  Ref  1,  p  135-140,  where  the  equation 

=  4nk6  /  I(Ni  +  N2)N3N4  -  (N3  +  N4)NiN2] 

•  6(Oi  +  O2  -  a3  -  04)6(ki  +  k2  -  k3  -  k4)dk2dk3dk4  ,  (25) 

=  N(k^),  is  found  to  be  valid  in  the  neighborhood  of  the  spectral  peak  at  k 
=  k^  (Ref  19).  See  also  Fig  4  of  Ref  9,  where  the  "overshoot  effect"  in  the 
spectral  density  development  is  shown. 

Recent  results  on  prediction  models  based  on  the  radiative  transfer 
equation  are  discussed  in  Ref  9,  p  721-722.  Barnett  (Ref  20,  21)  used 
finite-difference  methods  to  solve  the  radiative  transfer  equation  over  the 
entire  North  Atlantic  under  the  following  assumptions: 

3F  ♦  8F 

It  ^  *i  Ir  =  ^  ±  Ss  -  Se  ,  (26) 

i 

where  is  a  version  of  Phillips'  resonance  theory  modified  by  the 
observations  of  Section  5.1  of  Ref  9;  82  is  an  exponential  growth  term  based 
entirely  on  the  observations  described  in  Section  5.1  of  Ref  9;  85  represents 
computer-simplified  parameterizations  of  the  wave-wave  interactions;  and  83  is 
a  representation  of  wave  breaking  (limiting  growth)  based  on  the  equilibrium 
ideas  of  Phillips  as  described  in  Section  5.5  of  Ref  9.  No  assumption  of  a 

19.  M.S.  Longuet-Higgins ,  "On  the  nonlinear  transfer  of  energy  in  the  peak  of 
a  gravity-wave  spectrum:  a  simplified  model,"  Proc.  Roy.  Soc.,  A,  437,  1976, 
311-328. 

20.  T.P.  Barnett,  "On  the  generation,  description,  and  prediction  of  ocean 
wind  waves,"  Ph.D.  dissertation.  University  of  California,  San  Diego,  1966. 

21.  T.P.  Barnett,  "On  the  generation,  description,  and  prediction  of  ocean 
wind  waves,"  J.  Geophys.  Res.,  73,  1968,  513-529. 


fully  developed  spectrum  was  made.  The  same  approach  has  been  followed  by 
Ewing  (Ref  22) .  His  formulation  utilized  a  more  accurate  parametric  version 
due  to  Cartwright  of  the  term  85  and  a  higher-order  finite-difference 
approximation  for  the  advective  term  of  Eq  (19).  Ewing's  comparisons  were 
made  against  the  two-dimensional  spectrum,  the  first  such  verifications 
atteiapted.  The  results  gave,  in  Ewing's  words,  "adequate  estimates  of  the 
significant  height  and  one -dimensional  wave  spectrum.  The  standard  deviation 
of  all  the  computed  estimates  of  [the  significant  wave  (the  mean  of  the 

highest  third  of  all  waves  present)]  conqpared  to  measurements  is  about  0.6  m. 
Reliable  estimates  of  the  two-dimensional  wave  spectrum  were  only  achieved  in 
a  limited  region  at  the  high-frequency  end  of  the  spectrum."  This  later 
discrepancy  was  attributed  to  inadequate  specification  of  the  wind  field. 


22.  J.A.  Ewing,  "A  nuawrical  wave  prediction  for  the  North  Atlantic  Ocean," 
Deutsche  Hydrog.  2. .  24,  1971,  241-261. 
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III.  STATE  OF  THE  ART  IN  PRACTICAL  OCEAN  SURFACE  MODELING 


A.  INTRODUCTION 

An  important  empirical  relation  is  the  "mean  JONSWAP  spectrum,"  Eq  (27). 
This  is  discussed  in  Section  III.B  along  with  some  other  useful  empirical 
relations.  Using  Eq  (27)  allows  the  wind  speed  magnitude  to  be  considered; 
the  directional  behavior  of  the  wind  is  addressed  in  Section  II. E,  especially 
Eq  (17).  Surface  waves  in  a  variable  current  are  discussed  in  Section  III.C, 
where  direction  and  energy  density  relations  are  given  for  four  situations: 

(1)  upwelling,  (2)  lateral  spreading,  (3)  frequency- dependent  attenuation,  and 
(4)  refraction  of  surface  waves  by  a  current  shear.  Section  III.D  treats  the 
interaction  of  swell  and  short  gravity  waves.  Equation  (73)  represents  a  key 
result  for  imaging  by  a  remote  radar  sensor.  The  interaction  of  surface  and 
internal  waves.  Section  III.E,  is  our  best-understood  phenomenon.  Particular 
emphasis  is  given  to  the  approaches  of  Phillips  (Ref  1,  p  78-81,  and  Ref  23), 
Ko  (Ref  24),  and  Hughes  (Ref  17).  Finally  wave  attenuation  is  addressed  in 
Section  III.F. 

Weak  wave-wave  interactions  will  not  be  discussed  further  since  they  are 
covered  thoroughly  in  Ref  1,  p  27-32,  81,  135-140,  and  Ref  9,  p  682-684.  On 
the  other  hand,  long  wave-short  wave  interaction  is  still  a  controversial 
subject.  Garrett  and  Smith  (Ref  25)  update  the  survey  material  contained  in 
Ref  9,  p  690-691,  and  conclude  that  the  long  waves  can  grow  if  short  wave 


23.  O.M.  Phillips,  "On  the  interactions  between  internal  and  surface  waves," 
Phys.  Atmos,  and  Oceans,  9,  1973,  954-961. 

24.  J.E.  Lewis,  B.M.  Lake,  and  D.R.S.  Ko,  "On  the  interaction  of  internal 
waves  and  surface  gravity  waves,"  J.  Fluid  Mech. ,  63(4),  1974,  773-800. 

25.  C.  Garrett  and  J.  Smith,  "On  the  interaction  between  long  and  short 
surface  waves,"  J.  Phys.  Oceanogr..  6,  1976,  925-930. 


generation  is  correlated  with  the  long  wave  orbital  velocity.  Reference  26  is 
very  instructive  and  contains  some  tools  from  which  a  model  can  be  built. 
References  27  and  28  give  some  experimental  results;  the  latter  paper  uses 
swells  as  the  long  waves.  Valenzuela  and  Wright  (Ref  29)  use  the  radiative 
transfer  equation,  developed  for  gravity-capillary  waves,  but  their  work  is 
too  theoretical  for  easy  model  development. 


B.  THE  MEAN  JONSWAP  SPECTRUM 

Phillips  (Ref  1,  p  139-140)  and  Barnett  and  Kenyon  (Ref  9,  p  715-717) 
reviewed  the  "mean  JONSWAP  spectrum."  Fox  (Ref  30)  found  for  the  mean  JONSWAP 
spectrum  an  empirical  fit  to  results  from  the  Joint  North  Sea  Wave  Project 
(Ref  31)  as  follows: 


♦i(n)  =  yg*n-®  exp 


lWh)‘  ■ 


26.  O.M.  Phillips,  "The  dispersion  of  short  wavelets  in  the  presence  of  a 
dominant  long  wave,"  J.  Fluid  Mech.,  107,  1981,  465-485. 

27.  O.H.  Shemdin,  "Modulation  of  centimetric  waves  by  long  gravity  waves: 
Progress  report  on  field  and  laboratory  results,"  NATO  Conference  on  Turbulent 
Fluxes  Through  the  Sea  Surface,  Wave  Dynamics  and  Prediction,  A.  Farve  and  K. 
Hasselmann,  eds.,  NATO  Conference  Series  V,  vol.  1,  1977. 

28.  E.D.  Graves  and  O.H.  Shemdin,  "An  investigation  of  the  modulation  of 
capillary  and  short  gravity  waves  in  the  open  ocean,"  J.  Geophys.  Res., 

85(C9),  1980,  5019-5024. 

29.  G.R.  Valenzuela  and  J.W.  Wright,  "Modulation  of  short  gravity- capillary 
waves  by  longer-scale  periodic  flows — a  higher-order  theory,"  Radio  Science, 
14(6),  1979,  1099-1110. 

30.  M.J.H.  Fox,  "On  the  nonlinear  transfer  of  energy  in  the  peak  of  a 
gravity-wave  spectrum — I,"  Proc.  Roy.  Soc.,  A  348,  1976,  467-483. 

31.  K.  Hasselmann,  et  al,  "Measurements  of  wind  wave  growth  and  swell  decay 
during  the  Joint  North  Sea  Wave  Project  (JONSWAP),"  Herausgegeben  von  Deutsche 
Hydrograph.  Institut .  Reihe  A,  12,  1973. 
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m  =  D  for  n  <  n  and  m  =  m.  for  n  >  n  .  This  form  contains  five  paraowters: 
Y,  Yg*  frequency  of  the  spectral  peak  n^.  Note  that  Eq  (27)  is 

an  extension  of  Eq  (14)  to  all  frequencies.  For  the  mean  JONSWAP  spectrum,  Y^ 
=  3.3,  =  0.07  and  nij^  =  0.9.  Figure  17  of  Ref  9  gives  Y  «*  * 

function  of  the  nondimens ional  fetch 


X  =  d^g/U2, 


where  d^  is  the  fetch  and  the  U  is  the  wind  speed  in  the  forms 


Y  =  y(*)  =  0.076  X 

(see  also  Table  4.1  of  Ref  1)  and 

n  =  n  (x)  =  2n*3.5  x  . 

o  o 


In  Ref  1,  p  148,  161,  the  following  useful  relations  are  found: 


n  ~  2.2(g3/u?d  )^ 


*  /2  \  * 

u(2)  =  IT  (ir)‘  2r^^V>  ' 

o  \  o/  o 


where 


09  *2k  zt, 

f(k^2)  =  (k^2)*  J  t  ^  e  °  dt  , 


n^  =  1.1-10-2  u2/g  ^ 


K  -  0.42  is  Karman's  constant,  B  -  2*10  and  k  is  the  wavenumber  spectral 


C.  WAVES  ON  A  VARIABLE  CURRENT 


Waves  on  a  variable  current  are  treated  in  Ref  1,  p  74-78,  and  Ref  16. 
When  a  wave  train  encounters  a  current  in  which  the  surface  velocity  varies, 
the  excess  momentuin  flux  results  in  an  interchange  of  energy  between  waves  and 
current.  Suppose  that  the  current  strength  varies  along  the  direction  of  flow 
so  that  U  =  U(x)  in  Eq  (2);  a  situation  like  this  may  be  encountered  at  an 
estuarine  outflow  with  waves  incident  from  the  open  sea.  If  a  represents  the 
angle  between  the  local  wavenumber  and  the  current  in  the  x-direction,  then  Eq 
(1)  and  (2)  for  a  steady  wave  train  in  deep  water  give 

n  =  o  +  kU  cos  a  =  const.  =  =  8^%  ~ 

where  is  the  wave  intrinsic  frequency  in  the  open  sea,  where  U  =  0.  Also 
Eq  (3)  gives 

k  sin  Of  =  k^  sin  =  const.  (37) 

The  wavenumber  component  in  the  y-direction  is  conserved.  From  these  two 
expressions 

(gk)^  +  U[k*  -  k*  sin  ]^  =  a  (38) 

o  o  o 

With  a  current  in  the  same  direction  as  that  in  which  the  waves  are 
traveling,  U  >  0  and  the  wavenumber  k  is  reduced,  the  wavelength  increases, 
and  the  wave  paths  are  turned  away  from  the  current.  In  an  adverse  current, 
on  the  other  hand,  U  <  0,  and  the  waves  turn  into  the  current  as  they  shorten. 
Their  group  velocity  decreases  and  they  become  less  able  to  propagate  against 
the  current.  If  the  strength  of  the  latter  is  sufficient,  an  interesting 
phenomenon  of  blockage  occurs.  The  algebra  is  simplest  when  a  =  0  and  the 
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waves  are  running  directly  into  the  current.  Equation  (38)  then  becomes  a 
sioiple  quadratic  equation  in  k^,  or  equivalently  in  the  local  phase  velocity  c 
=  (g/k)^»  with  the  solution 


■■•"('•S' 


which  has  a  critical  point  at  U  =  -%c  =  where  energy  can  no  longer  be 

propagated  against  the  current.  This  is  a  kinematical  limit;  the  waves  break 
before  this  point.  In  summary,  for  waves  in  the  presence  of  an  adverse 
current,  the  waves  turn  toward  the  stream.  Their  local  wavelength  decreases 
and  the  energy  density  increases  as  the  waves  are  pushed  toward  saturation  and 
possible  blockage.  In  a  favorable  current  they  turn  away  from  the  current, 
the  local  wavelength  increases,  the  energy  density  is  reduced,  and  the  waves 
become  less  saturated. 

Equation  (5)  gives 


E(U  +  ^c)c  =  const.  =  VE  c* 

o  o 


which  was  originally  derived  by  Longuet-Higgins  and  Stewart  (Ref  7)  from  Eq 
(7).  In  the  absence  of  wave  breaking,  the  local  amplitude  is  given  by  (see 
Fig  3.6  of  Ref  1) 


Ic(c  +  2U)] 


=(U' 


When  the  convergence  of  the  current  is  balanced,  not  by  upwelling  as 
above,  but  by  lateral  spreading. 


+  §Y  ,  0 

8x  8y  “  • 


If  the  flow  is  symmetrical  about  the  line  y  =  const. ,  then  on  this  line 


(43) 


Now  assume  Eq  (16)  applies.  Suppose  that  the  wind-generated  waves  are 
unidirectional  and  move  in  direct  opposition  to  the  local  current.  In  a 

frequency  band  do  in  still  water,  the  wave  energy  is  initially,  via  Eq  (16), 

♦^(a)do  =  Yg^o~®do  .  (44) 

At  the  point  where  the  adverse  current  is  the  greatest,  (U  =  -U  ) ,  the 

m 

frequency  of  this  band  is  o'  and  the  wave  energy  is 

♦i(o)  =  Y«^(ff’)"®  .  (45) 

As  the  current  decreases  to  zero,  the  frequency  of  the  band  returns  to  o  and 

the  wave  energy  is  reduced  by  the  expansion  to  In  this  part  of  the 

motion,  no  energy  is  lost  by  wave  breaking.  Equation  (40)  gives 


(46) 

(47) 

(48) 

(49) 

(50) 


A  slowly  varying  current  case  is  treated  in  Ref  32. 


Refraction  of  surface  waves  by  shear  is  probably  encountered  more 
frequently  in  the  ocean  than  the  situation  just  described,  particularly  at 
oceanic  fronts,  where  a  contrast  in  surface  temperature  is  necessarily 
accompanied  by  shear.  Since  surface  waves  are  dispersive,  refraction  effects 
by  currents  are  stronger  than  they  are  in  nondispersive  waves. 

The  wave  patterns  can  be  analyzed  in  a  manner  similar  to  that  given 
above.  If  a  represents  the  angle  between  the  local  wavenumber  and  the  current 
and  U  =  U(y)  represents  the  shearing  surface  current,  Eq  (1)  and  (2)  give 

a  +  kU  cos  a  =  (51) 

but  now  k  =  k(y)  in  Eq  (3)  gives 

k  cos  a  =  k^  cos  =  const.  (52) 

From  these  expressions  and  the  dispersion  relation  o  =  (gk)^,  it  follows  that 


cos  a 


gk  cos  a 
_  o _ o 

(a  -  k  U  cos  a  )^  ’ 

o  o  o 


(53) 


which  is  necessarily  less  than  or  equal  to  unity.  If  cos  <  0  and  the  waves 
approaching  the  current  are  opposed  to  it,  they  will  be  turned  toward  the 
normal.  If  cos  >  0,  however,  they  are  turned  away,  the  effect  being 
augmented  beyond  simple  convection  by  the  dispersive  nature  of  the  waves. 

Just  as  the  local  wavelength  increases  as  a  result  of  the  stretching,  so  does 
the  wave  velocity.  The  maximum  penetration  of  the  waves  into  the  current 
system  is  reached  when  a  =  0  and 


32.  D.  Holliday,  "Nonlinear  gravity-capillary  surface  waves  in  a  slowly 
varying  current,"  J.  Fluid  Mech..  ^(4),  1973,  797-802. 


(54) 


y  1  -  (cos 

c  cos  a 

o  o 

The  waves  are  then  refracted  back  by  the  stream.  When  the  initial  angle  of 

incidence  between  the  waves  and  the  current,  a^,  is  relatively  small,  the 

stream  current  needed  to  turn  back  the  waves  is  also  small,  approximately 

a^c  /4,  where  a  is  measured  in  radians, 
o  o  o 

The  distribution  of  wave  energy  density  can  be  found  from  Eq  (5).  Since 

U  =  U(y),  the  flux  of  wave  action  normal  to  the  current  is  independent  of  y: 

(Ec  /a)  sin  a  =  const.  (55) 

E 

or,  since  c^  =  g/2a, 

E  sin  a  «  =  (a^  -  U(y)k^  cos  of^I^  (56) 

from  Eq  (51)  and  (52),  with  again  a  theoretical  singularity  in  E  at  or  =  0  when 
Eq  (54)  is  satisfied.  Note  that  from  Eq  (53)  and  (56) 

E  sin  2or  =  const.  (57) 

Waves  refracted  toward  the  normal  to  the  current  are  reduced  in  energy 
density;  in  those  refracted  away  from  the  normal,  the  local  energy  density  is 
augmented.  The  same  is  true  for  the  mean  square  slope,  which  is  proportional 
to  Ek^ : 


Ek^ 


sin  2a  ^  sin  4o 


(58) 


“o  o  sin  2a  +  ^  sin  4a 
so  that  wave  trains  refracted  back  from  the  current  are  driven  toward  local 
saturation  as  the  turning  line  a  =  0  is  approached. 

How  may  these  variations  be  reflected  in  the  SAR  (synthetic  aperture 
radar )^  imagery,  in  which  a  given  wavenumber  on  the  surface  is  sampled, 


^See  appendix. 
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not  a  given  wave  train?  If  the  instruunent  is  looking  directly  across  the 
current,  a  =  =  n/l\  the  wavenumber  and  energy  density  are  unchanged  and  there 

is  little  reason  to  expect  that  the  stream  will  be  detected  at  all.  If  the 
current  U  has  a  positive  component  along  the  direction  of  observation,  as  U 
increases  with  range,  for  example,  it  will  sample  waves  that  are  closer  to 
local  saturation;  but  beyond  the  line  of  maximum  wave  penetration,  the  return 
signal  would  be  expected  to  be  much  less.  If  U  has  a  negative  component  along 
the  direction  of  observation,  outwardly  travelling  waves  will  be  refracted 
toward  the  normal;  their  mean  square  slope  decreases  and  they  become  less 
saturated,  giving  a  reduced  signal  intensity,  provided  the  wind  regeneration 
distance  is  large  enough.  The  behavior  of  approaching  wave  trains  can  be 
inferred  similarly. 

D.  INTERACTION  OF  SWELL  AND  SHORT  GRAVITY  WAVES 

We  turn  to  Phillips  (Ref  16)  for  a  discussion  of  the  modification  of 
short  waves  by  swell.  The  short  wave  components  are  converted  and  distorted 
by  longer  waves  or  swell.  Short  wave  energy  and  amplitude  are  in  fact 
concentrated  near  the  crests  of  the  long  waves  and  somewhat  depleted  in  the 
troughs.  In  an  actively  generated  wind-wave  field,  the  short  components  may 
then  be  saturated  only  locally,  near  the  long  wave  crests,  and  somewhat 
undersaturated  elsewhere.  The  long  waves  propagate  faster  than  the  short 
waves  so  that  a  given  group  of  short  waves,  being  overtaken  by  a  swell  crest, 
may  become  saturated  with  local  sporadic  breaking,  then--as  the  crest  moves 
past--become  unsaturated  as  it  slips  into  the  trough.  It  continues,  however, 
to  receive  energy  from  the  wind;  and,  if  the  input  is  sufficient,  it  is  ready 


again  for  more  occasional  wave  breaking  just  before  and  at  the  next  crest  of 
the  swell.  This  local  saturation  at  the  crests  of  the  swell  thus  produces  a 
modulation  in  the  return  signal  of  a  SAR  that  enables  the  long  wave  pattern  to 
be  discerned  if  its  wavelength  is  greater  than  about  twice  the  discrimination 
distance.  This  effect  is  enhanced  somewhat  by  a  stationary  wave  effect  when 
short  waves,  whose  wavelength  varies  with  respect  to  the  phase  of  the  swell, 
are  viewed  with  the  fixed  radar  wavelength  signal.  The  steepness  of  the 
spectrum  of  Eq  (15)  provides  further  enhancement.  Field  studies  by  Evans  and 
Shemdin  (Ref  28)  have  demonstrated  this  modulation  in  radar  return  signals  and 
associated  them  clearly  with  variations  in  short  wave  structure  produced  by 
longer  waves . 

The  modulation  of  short  wave  structure  by  swell  was  first  analyzed 
correctly  by  Longuet-Higgins  and  Stewart  (Ref  6),  but  their  analysis  was 
limited  to  conditions  in  which  the  swell  amplitude  is  less  than  the  short  wave 
wavelength.  Numerical  calculations  that  are  free  of  this  restriction  were 
performed  by  Longuet-Higgins  (Ref  33),  and  if  (as  is  the  situation  of  interest 
here)  the  ratio  of  long  to  short  wavelength  is  large,  simple  asymptotic 
methods  can  be  used.  The  details  of  this  are  given  by  Phillips  (Ref  26). 

These  analyses  were  concerned  with  individual  wave  trains.  For  the 
present  purpose  they  must  be  extended:  (a)  to  allow  for  an  arbitrary 
direction  of  propagation  of  the  short  waves  relative  to  the  swell,  (b)  to 
express  the  results  of  the  interaction  in  terms  of  the  local  short  wave 
spectral  density,  and  (c)  to  determine  the  magnitude  of  the  modulations  in 
short  wave  spectral  density  of  a  fixed  wave  number  on  the  surface,  not  on  a 
given  energy  path  along  which  the  wavenumber  varies. 


33.  M.S.  Longuet-Higgins,  "The  instabilities  of  gravity  waves  of  finite 
amplitude  in  deep  water.  Part  I:  Superharmonics,"  Proc.  Roy.  Soc.  A  360, 
1978,  471-488. 


21 


The  geometry  is  illustrated  in  Fig  1.  Axes  x  and  y  are  chosen  in  the 
direction  of  swell  propagation  and  along  the  crests,  respectively;  (|i  is  the 
angle  between  the  direction  of  observation  and  the  direction  of  swell  propa¬ 
gation,  and  0^  is  the  direction  of  the  local  wind  relative  to  the  viewing 
direction.  In  a  frame  of  reference  moving  with  the  swell,  let  its  surface 

displacement  be  ^  =  A  cos  Kx  =  A  cos  x*  The  short,  locally  generated  wind 

s 

waves  are  riding  over  the  swell;  the  local  wavelength  of  a  wave  train  and  its 
direction  of  propagation  vary  as  a  result  of  the  variable  orbital  velocities 
and  vertical  accelerations  produced  by  the  swell.  In  this  frame  of  reference 
the  near-surface  velocity  is 

U  =  -C(l  -  AK  cos  X)  (59) 

in  the  positive  x-direction,  where  C  and  AK  represent  the  phase  speed  and 
slope  of  the  swell,  respectively.  The  effective  gravitational  acceleration 
experienced  by  the  short  waves  is 

g'  =  g(l  -  AK  cos  x)  (60) 

as  shown  in  Ref  26. 

LOCAL  WIND 


Figure  1 .  Short  wave  trajectories  riding  on  a  swell. 


to  the  lowest  order  in  AK.  The  modulations  in  wavenumber  are  strongest  when  ^ 
=  0  and  disappear  when  cos  (|>  =  c/2C  «  1,  that  is,  when  the  short  waves  are 
running  almost  along  the  crests  of  the  swell.  It  is  interesting  to  note  that 
short  waves  propagating  exacting  along  the  swell  crests  are  deformed  somewhat 
since  the  effective  gravitational  acceleration  experienced  by  the  short  waves 
is  smaller  there  and  the  propagation  speed  is  smaller  than  in  the  troughs. 
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Between  long  wave  crests,  when  the  short  waves  are  unsaturated,  the 
dynamics  of  the  interaction  are  expressed  most  simply  by  Eq  (18).  For  short 
gravity  wavelengths  of  the  order  of  30  cm,  energy  input  from  the  wind  from  one 
swell  crest  to  the  next  can  be  ignored  if,  as  is  generally  the  case,  of  Eq 
(67)  is  substantially  larger  than  the  swell  period;  is  the  time  scale  over 
which  the  energy  density  decreases  by  a  factor  e  in  the  absence  of  dissipative 
processes  such  as  wave  breaking  and  when  c  >  10  u^  or  so,  and  is  given  by 

~  20c2/ou|  (67) 

or 

~  3(c^/u^)^  (wave  periods)  .  (67) 

In  a  frame  of  reference  moving  with  the  swell,  the  distribution  of  N(k)  is 
steady  so  that 

(c  +  y)‘VN(k)  =  0  .  (68) 

8 

Consequently  N(k)  is  constant  along  ray  trajectories,  though  k  =  k(x),  and  so 
if  the  local  wind  field  is  uniform,  it  is  constant  throughout.  Thus,  in 
particular, 


dN  ^  3N  ,  3N 

dx  3x  3k.  3x 


(69) 


where  3N/3x  expresses  the  variation  at  a  fixed  wavenumber  k  as  sampled  by  the 
SAR.  Furthermore,  since  3k2/3x  =  0, 
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in  terms  of  the  polar  coordinates  of  Fig  1.  The  amplitude  of  the  modulations 
in  short  wave  action  spectra  at  a  fixed  wavenumber  is  then 


from  Eq  (65)  and  is  proportional  to  the  swell  slope  AK. 

Further  progress  is  dependent  on  some  knowledge  of  the  N(k).  If,  near 
the  swell  crests,  the  locally  generated  wind  waves  are  close  to  saturation  at 
the  wavenumbers  providing  the  backscattering,  with  local  energy  input  from  the 
wind  being  balanced  largely  by  sporadic  small-scale  wave  breaking,  the 
spectral  density  of  the  short  waves  is  given  by  Eq  (15),  where  f(6)  is 
symmetrical  about  the  wind  direction.  The  corresponding  action  spectral 
density  is 

N(k)  =  i|»(k)/o  =  Df(0)g‘V’/^  ,  (72) 

so  that  from  Eq  (65),  (70),  and  (71),  the  amplitude  of  the  modulations  6i|i  in 
short  wave  spectral  density  at  a  fixed  wavenumber  referred  to  the  mean  is 

^  ^  ^ 

«1>  ■  N  "  F 

=  [|  cos  ((|»  -  0^)  -  II  sin  ((>  -  0^)]g(C,c,4>)(AK).  (73) 

According  to  Bragg  scattering  theory,  the  amplitude  of  the  modulations  in 
return  signal  intensity  relative  to  the  mean  is  equal  to  6()</tli  and  so,  from 
this  expression,  is  proportional  to  the  swell  slope  AK. 
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The  result  given  above  shows  how  the  swell- induced  modulations  depend  on 
the  viewing  direction  ^  relative  to  the  swell  and  on  the  distribution  of 
short,  locally  generated  waves  about  the  wind  direction  6^.  The  angle 
between  the  direction  of  propagation  of  the  swell  and  the  direction  of  viewing 
can  be  determined  directly  from  the  SAR  imagery,  but  the  direction  of  the 
local  wind  cannot.  This  requires  either  a  simultaneous  scatterometer 
BMasurement  at  shorter  wavelengths  or  some  independent  information  about  the 
local  meteorology.  The  phase  speed  C  of  the  swell  can  be  deduced  from  the 
wavelength  of  the  observed  swell  pattern  and  the  mean  phase  speed  c  of  the 
scatterers  by  a  knowledge  of  the  Bragg  wavelength.  Equation  (73)  does  not 
explicitly  involve  the  local  wind  speed,  though  clearly  there  are  limits  to 
the  range  of  wind  speeds  over  which  it  is  applicable.  The  wind  must  be 
sufficiently  strong  to  generate  locally  short  waves  at  the  Bragg  scattering 
wavelength  but  not  so  large  that  the  dominant  locally  generated  waves  have 
orbital  velocities  comparable  with  those  of  the  swell.  In  this  latter  case, 
the  modulation  of  the  short  waves  by  the  longer,  locally  generated  waves  may 
be  dominant  but,  being  on  a  scale  smaller  than  the  resolution  distance  of  the 
SAR,  will  not  be  observed.  Beal  (Ref  34)  gives  a  range  of  2  m/s  <  < 

10  m/s  (where  is  the  wind  speed  at  10  m  above  the  mean  sea  surface)  over 
which  swell-induced  modulations  have  been  detected. 

An  ideal  observing  situation  would  be  one  in  which  0^  =  ~  0®  with  swell 

moving  in  the  viewing  direction  and  a  following  (or  opposing)  local  wind.  The 
relative  modulation  amplitude  then  reduces  simply  to 

^  =  4.5AK  .  (74) 


34.  R.C.  Beal,  "Spaceborne  imaging  radar:  monitoring  of  ocean  waves," 
Science.  208,  1980,  1373-1375. 
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As  ^  increases  fron  zero,  since  C  »  c  in  the  last  factor  of  Eq  (73),  the 
Modulation  amplitude  decreases,  reaching  a  minimum  near  n/2.  This  is 
consistent  with  Beal's  remark  (Ref  34)  that  the  swell  can  be  imaged  success¬ 
fully  when  there  is  a  substantial  component  of  the  swell  traveling  along  the 
line  of  sight  of  the  radar. 

The  first  factor  in  Eq  (73)  involves  the  spectral  distribution  of  the 
short  locally  generated  wind  waves  that  provide  the  surface  scatterers,  the 
first  term  representing  the  consequences  of  variations  in  wavenumber  magnitude 
and  the  second  in  local  short  wave  direction.  The  directional  distribution  of 
the  short  gravity  components  in  a  wind-generated  wave  field  is  discussed  in 
Section  II. E.  Unless  the  directional  distribution  of  these  waves  is  very 
narrow,  however,  the  second  term  is  likely  to  be  numerically  rather  smaller 
than  the  first,  especially  when  (^  -  0^)  is  small.  Under  these  circumstances, 
a  simplified  form  of  Eq  (73)  may  be  useful  in  estimating  the  swell  slope  AK 
directly  from  satellite  measurements. 

In  an  experimental  paper  Evans  and  Shemdin  (Ref  35)  observe  that  the 
modulation  of  short  waves  by  swell  is  found  to  be  strong  enough  to  be  an 
iaiportant  component  of  SAR  image  formation  for  the  swell.  Phillips  (Ref  1, 
p  97-98)  presents  some  further  analysis  of  the  action  of  swell  moving  across 
the  sea  surface.  He  observes  that  if  the  augmented  drift  U  of  Eq  (59)  near 
the  swell  crests  approaches  the  local  phase  speed  c  of  the  short  waves,  the 
latter  may  be  suppressed  entirely.  This  effect  was  apparently  observed  first 


35.  D.D.  Evans  and  D.H.  Shemdin,  "An  investigation  of  the  modulation  of 
capillary  and  short  gravity  waves  in  the  open  ocean,"  J.  Geophys.  Res.,  85, 
1980,  5019-5024. 


by  Hitsuyasu  (Ref  36)  and  analyzed  by  Phillips  and  Banner  (Ref  37),  who  also 
conducted  further  laboratory  neasureaents . 

L 

The  dispersion  relation  for  short  waves  is  a  =  (g'k)  [see  Eq  (4)  and 
(60)],  and  for  short  waves  in  the  frasie  of  reference  of  the  swell,  Eq  (1)  and 
(2)  give 

a  +  k(U  -  C)  =  const.  (75) 

If  c^  is  the  short  wave  speed  at  the  point  X  =  n/2,  where  the  surface 
displacement  of  the  swell  is  zero  (remember  that  x  =  0  at  the  swell  crest)  and 
k^  is  the  corresponding  wavenumber, 

(g'k)^  +  k(U  -  C)  =  (gk^)^  -  k^C 

=  k^(c^  -  C) 
o  o 

=  eCCq  -  C)/cJ  .  (76) 

This  is  a  quadratic  in  k  or  equivalently  in  c  =  (g'/k)’,  and  since  U  = 
and  cos  x  =  CAK  cos  Xi  it  can  be  written  as 

im--  'ak°c<,s  x)(§)'  * 

where  a  =  c^/C  <<  1.  Since  c/C  a  when  x  n/2,  the  appropriate  root  is 
c/C  =  a(l  -  Ak  cos  x)  •  (78) 

36.  H.  Mitsuyasu,  "Interactions  between  water  waves  and  wind  (I)," 

Rep.  Res.  Appl.  Mech.,  Kyushu  University,  14,  1966,  67-88. 

37.  O.M.  Phillips  and  M.L.  Banner,  "Wave  breaking  in  the  presence  of  wind 
drift  and  swell,"  J.  Fluid  Mech..  66,  1974,  625-640. 


Note  that  the  phase  speed  c  of  the  short  waves  is  reduced  as  the  crest  is 
approached,  while  the  aabient  drift  q  increases.  Consequently  the  ratio  q/c 
increases  and  the  Baxiauai  height  of  an  unbroken  wavelet  decreases.  The 
■axiauai  asqtlitude  that  wavelets  can  have  when  they  are  at  the  point  of 
incipient  breaking  at  the  swell  crest  is,  from  Eq  (3.9.6)  of  Ref  1, 


c 


■ax 


(2g*)''(c^  -  q^)2  , 


(79) 


where  c^  =  Cq(1  ~  AK)  is  the  phase  speed  of  the  short  wave  at  the  swell  crest, 
and  q^  is  the  augnented  value  of  the  drift  at  this  point,  given  in  Eq  (3.9.4) 
of  Ref  1  as 


q^  =  (C  -  U^)  -  I(C  -  U^)2  -  q^(2C  -  q^)]^  . 


(80) 


The  ratio  r  of  the  maximum  wavelet  amplitude  to  its  maximum  in  the  absence  of 
the  swell  is,  from  Eq  (3.9.6)  of  Ref  1  and  Eq  (79), 


r 


g' 


(81) 


which  is  an  algebraic  function  of  the  three  parameters  a  =  c^/C,  AK  =  U^/C  and 
y  =  Qq/C.  Measurements  by  Mitsuyasu  (Ref  36)  and  Phillips  and  Banner  (Ref  37) 
on  the  suppression  of  short  wind-generated  waves  in  the  laboratory  by  longer, 
mechanically  generated  swell  gave  good  agreement  with  Eq  (81)  over  the  rather 
limited  range  of  conditions  studied. 
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E.  INTERACTION  OF  SURFACE  AND  INTERNAL  WAVES 


1 .  Introduction 

Qualitative  relations  for  the  behavior  of  short  surface  waves  interacting 
with  internal  waves  are  found  in  Ref  1,  p  78-81,  230-232,  Ref  23,  and  Ref  38. 
Lewis,  Lake,  and  Ko  (Ref  24)  do  a  one-dinensional  theoretical  analysis  of  a 
tank  experiment.  They  demonstrate  that  the  normalized  amplitude  and  slope 
iBodulations  increase  with  the  interaction  distance  and  that  the  maximum 
interaction  effect  occurs  when  the  phase  speed  of  the  internal  wave  equals  the 
group  speed  of  the  surface  wave.  Hughes  (Ref  17)  produced  a  theoretical  model 
that  showed  promising  agreement  with  the  experimental  results  of  Ref  39. 
Holliday  (Ref  32)  shows  that  the  "wave  barrier"  noted  by  Gargett  and  Hughes 
(Ref  38)  for  infinitesimal  gravity  waves  on  a  slowly  varying  current  from  an 
internal  wave  is  removed  by  considering  finite  amplitude  effects. 

This  section  considers  in  detail  the  results  of  Phillips  (Ref  1,  p  78-81, 
and  Ref  23),  Lewis,  Lake,  and  Ko  (Ref  24);  and  Hughes  (Ref  17).  First,  the 
Phillips  approach  is  given.  It  considers  the  blockage  of  surface  waves  by 
interaction  with  internal  waves.  In  the  absence  of  blockage,  when  energy 
packets  move  continuously  across  the  surface,  he  finds  that  a  steady  state  can 
be  obtained  in  which  the  energy  ratio  for  surface  is  given  by  Eq  (88).  He 
then  proceeds  to  show  that  when  the  wide  spectral  range  of  wind-generated 
waves  is  considered,  the  blockage  singularity  at  U  +  ^c  -  C  =  0  in  Eq  (88) 
vanishes  and  Eq  (93)  results.  Finally,  he  considers  the  problem  from  the 


38.  A.E.  Gargett  and  B.A.  Hughes,  "On  the  interaction  of  surface  and  internal 
waves,"  J.  Fluid  Mech.,  ^(1),  1972,  179-191. 

39.  B.A.  Hughes  and  H.L.  Grant,  "The  effect  of  internal  waves  on  surface  wind 
waves.  1:  Experimental  measurements,"  J.  Geophys.  Res.,  83.  1978,  443-454. 


point  of  view  of  weak  wave-wave  interactions  and  finds  the  resonance  condition 
of  Eq  (96). 


Kext  the  Ko  approach  is  reviewed.  To  nodel  the  Lewis,  Lake  and  Ko  (Ref 
24)  tank  experisient  he  begins  with  the  appropriate  one-dinensional  foms  of  Eq 
(1)  and  (7),  applies  the  characteristic  Eq  (102)  and,  through  a  perturbation 
analysis,  obtains  Eq  (129)  and  (130).  Both  the  theoretical  and  experimental 
results  for  slope  and  amplitude  modulations  show  that,  when  normalized  by  the 
ratio  of  induced  surface  current  to  internal  wave  phase  speed  (ie,  u^/C),  the 
resonant  case  (ie,  internal  wave  phase  speed  C  matches  the  surface  wave  group 
speed  Cgg)  Bodulation  magnitudes  are  dependent  only  on  the  interaction 
distance  and  are  independent  of  the  magnitudes  of  the  other  physical 
parasieters.  Likewise,  the  maximum  normalized  slope  and  amplitude  interaction 
effects  are  found  to  occur  when  C  =  c 


go 

Finally  Hughes'  approach  is  treated.  The  Hughes  (Ref  17)  model  is 
designed  for  comparison  with  the  experimental  results  of  Ref  39.  He  begins 
with  the  appropriate  forms  of  Eq  (1),  (2),  and  (4)  with  a  source  term  added  to 
the  latter.  This  is  Eq  (133),  which  he  manipulates  to  a  form  of  Eq  (8)  in  Eq 
(134),  with  the  source  term  given  by  Eq  (135),  which  is  a  form  of  the  first 
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three  terms  1  S.  of  Eq  (22)-(24).  After  some  manipulation  he  arrives  at  Eq 
i=l  ^ 

(140),  which  is  a  simplified  form  of  the  radiative  transfer  equation  in  terms 
of  the  action  spectrum.  From  this,  to  the  first  order  in  |U|/C,  he  obtains  Eq 
(143c).  Expressing  the  internal  wave  surface  speed  U  as  a  sum  of  different 
internal  wave  wavenumber  components  [Eq  (144)],  he  then  obtains  Eq  (145)  for 
the  wave  displacement  spectrum  i|>(kj,k2),  from  which  he  obtains  some  quantities 
which  he  can  compare  with  the  experimental  results  of  Ref  39.  The  resulting 
comparison  of  theory  and  experimental  results  is  encouraging. 
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2.  Phillips'  Approach 


The  Phillips  approach  to  the  modulations  in  surface  waves  induced  by 
internal  waves  is  given  in  Ref  1  and  23  (p  78-81).  When  long  internal  waves 
of  the  kind  described  in  Section  5.3  of  Ref  1  propagate  with  velocity  C  along 
the  oceanic  thermocline,  the  orbital  motion  near  the  surface  is  oscillatory, 
in  essence  horizontal  and  of  the  form  U(K*(x  -  Ct)],  where  K  is  the  wavenumber 
of  the  internal  wave.  As  far  as  short  surface  waves  are  concerned,  provided 
that  their  wavenumber  k  >>  K,  the  disturbance  current  induced  by  the  internal 
wave  is  independent  of  depth  within  the  range  |kz|  <  1  and  has  a  horizontal 
length  scale  (the  internal  wave  wavelength)  that  is  very  much  greater  than 
their  own  wavelength.  The  methods  of  Section  3.6  of  Ref  1  can  then  be  applied 
to  such  a  train  of  short  surface  waves  as  they  are  modified  by  the  moving 
current  distribution  y(K*(x  -  Ct)),  supposed  known. 

The  algebra  is  simplest  when  the  two  wave  trains  are  collinear.  In  a 
frame  of  reference  moving  with  the  internal  wave,  the  pattern  of  surface  waves 
when  steady  gives,  from  Eq  (1)  and  (2), 


a  -  k(C  -  U)  =  const.  =  a  -  k  C  , 

o  o 


where  =  (gk^)^  is  the  surface  wave  intrinsic  frequency  at  the  phase  point 
of  the  internal  wave,  where  the  induced  current  vanishes.  This  equation  can 
be  written  as  a  quadratic  in  c/c^  with  the  solution 


^  =  I2(c^  -  C)]"^  c„  +  (c_  -  C) 


4U(c  -  C)' 
'  *  (c„  -  2C)5 


which  is  real,  provided 


If  at  some  point  equality  occurs, 


c 


and  it  can  be  shovm  that  ^c  +  U  -  C  =  0.  Relative  to  the  internal  wave 
pattern,  the  energy  flux  of  the  surface  wave  vanishes.  Note  that  if  the 
undisturbed  group  speed  approaches  the  phase  speed  C  of  the  internal  waves,  it 
appears  from  Eq  (84)  that  an  arbitrary  small  adverse  current  is  sufficient  to 
induce  blockage. 

The  distribution  of  wave  energy  can  be  found  from  the  action  conservation 
principle.  In  a  frame  of  reference  moving  with  the  speed  C  of  the  current 
pattern,  Eq  (5)  becomes 

-  S  [("  *  'g  -  ‘=>5]  =  0  •  <“> 

For  a  steady  wave  train,  since  a  =  g/c, 

Ec(U  +  ^c  -  C)  =  const.  =  E  c  (^c  -  C)  ,  (87) 

00  o 

where  E^  is  the  energy  density  at  the  phase  point  of  the  internal  wave,  where 
U  =  0.  As  U  varies  with  respect  to  the  phase  of  the  internal  wave,  the  local 
phase  speed  c  varies  also  as  specified  by  Eq  (83) ,  with  corresponding  modula¬ 
tions  in  the  energy  density.  Near  a  blockage  point,  U  +  ^c  -  C  ->  0  and  from 
Eq  (87),  the  steady-state  solution  requires  an  indefinitely  large  wave  energy 
density  and  evidently  fails.  From  a  given  init_al  state,  a  time-dependent 
solution  can,  of  course,  be  calculated  numerically  from  Eq  (1)  and  (86).  In 
the  absence  of  blockage,  when  energy  packets  move  continuously  &  ss  the 
surface,  a  steady  state  can  be  obtained  in  which  the  energy  ratio  of  the 
surface  waves  is 
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(88) 


F  c  (%c  -  C) 

a  _  o  o _ 

“  c(U  +  -  C) 


in  agreement  with  an  expression  given  by  Gargett  and  Hughes  (Ref  38)  for  this 

case.  The  singularities  in  the  energy  density  that  occur  in  Eq  (88)  if  the 

denominator  vanishes  are  clearly  associated  with  the  kinematic  limits, 

described  previously,  when  the  energy  propagation  speed  relative  to  the  moving 

pattern  vanishes  or  when  the  phase  speed  itself  vanishes  and  the  wavenumber 

tends  to  infinity.  Under  either  of  these  conditions,  the  steady-state 

solution  will  fail.  The  time-dependent  Eq  (86)  becomes,  at  the  points  on  the 

surface  where  U  +  c  -  C  =  0, 

8 


§E  _  au 

at  "ax' 


(89) 
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Accordingly,  E  =  E^e  ,  where  E^  is  the  energy  density  at  the  blockage  point 
at  time  t  =  0  and  a  represents  the  rate  of  surface  divergence  at  this  point. 

If  a  is  negative  (at  a  convergence  point)  the  energy  density  increases 
monotonically . 

The  analysis  to  this  point  has  been  concerned  with  a  single  train  of 
waves,  slowly  varying  in  space  and  time  as  a  result  of  the  interactions  with 
the  current  distribution.  Wind-generated  waves,  however,  extend  over  a  wide 
spectral  range,  and  the  previous  results  can  be  extended  to  describe  this 
situation.  When  the  surface  waves  are  sufficiently  two-dimensional,  simple 
analytical  results  can  be  derived  directly  from  the  solutions  given  above.  If 
F(k)  represents  the  spectral  density  in  the  wave  field,  then  the  energy 
contained  in  the  wavenumber  band  dk  is  F(k)dk.  The  total  energy  content  (per 
unit  area)  is  specified  in  the  steady  state  in  terms  of  the  energy  content 
F(k)dk  of  the  band  in  its  undisturbed  state  by  Eq  (88): 
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F(k)dk 


c  (%c  -  C) 

o  o 


F  (k)dk  c(U  +  -  C) 

o  o 

Note  that  the  width  of  the  band  varies  with  x  as  well  as  the  energy  density. 
The  variation  c;in  be  found  from  Eq  (82)  in  the  form 

kKg/k)**  +  U  -  C]  =  k^((g/k^)^  -  C]  .  (91) 


Taking  differential  increments  in  the  wavenumber  k  and  k^  for  a  given  U  and  C 
gives 

..  -  C 

^  = _ ?  _  (92) 

dk  -*  U  -  C  • 

o 

The  combination  of  Eq  (90)  and  (92)  leads  to  the  surprising  simple  result 


F(k)  _ 


F(k  ) 
o 


=  -  =  (k) 


An  interesting  property  of  this  expression  is  that  the  singularity  in  Eq  (90), 
when  the  net  group  speed  vanishes,  has  disappeared  from  Eq  (93).  The  reason 
for  this  is  that  the  range  of  wavenumbers  about  the  critical  value  at  the 
kinematic  limit  becomes  large,  as  does  the  energy  content,  so  that  the  energy 
density  in  wavenumber  space  remains  finite.  The  singularity  when  c  0 
remains,  though  for  situations  of  interest  in  this  section,  this  limitation  is 
not  encountered. 

It  is  instructive  to  examine  briefly  the  same  example  from  the  point  of 
view  of  the  resonant  interactions  discussed  in  Section  2.8  of  Ref  1.  Two 
surface  wave  trains  with  wavenumbers  ki,k2  can  form  a  resonant  triad  with  a 
col linear  internal  wave  with  wavenumber  ks  =  K,  provided  the  conditions 


ki  =  ka  +  kg  ,  Oi  =  02  +  03  , 
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where  O3  is  the  frequency  of  the  internal  wave,  are  satisfied.  Since  the 
wavelength  of  the  internal  wave  is  very  much  greater  than  those  of  the  surface 
waves,  k3  =  K  <<  ki,k2,  and  the  difference  kj  -  k2  =  6k  between  the  two 
surface  wavenumbers  is  very  small.  If  Oj  -  02  =  6a,  the  resonance  conditions 
of  Eq  (94)  are 

6k  =  k3  ,  6a  =  a3  (95) 


so  that 


6a  _  aa 
6k  k3 


or  c 

g 


=  C  . 


(96) 


For  resonance,  the  phase  speed  of  the  internal  wave  is  equal  to  the  group 
speed  of  the  surface  wave.  This  is  identical  to  the  condition  for  blockage 
with  an  internal  wave  of  small  amplitude.  Under  these  conditions  a 
steady-state  solution  is  not  possible.  Energy  resident  initially  in  a  uniform 
surface  wave  train  with  wavenumber  kj  flows  to  a  neighboring  wavenumber  kj  - 
6k  (and  to  the  internal  wave).  The  wave  train  develops  modulations  on  the 
scale  of  the  internal  wave,  their  amplitude  continually  increasing  with  time. 

The  resonant  generation  of  internal  waves  by  nonlinear  interactions  with 
surface  waves  has  been  studied  by  Brekhovskikh,  et  al  (Ref  40).  This  theory 
predicts  that  an  initially  uniform  surface  wave,  in  the  presence  of  an 

t 

internal  wave  such  that  Eq  (96)  is  satisfied,  will  develop  growing  modulations 
as  an  additional  Fourier  component  of  the  surface  wave  with  an  adjacent  wave- 
number  (the  third  member  of  the  triad)  is  generated.  The  two  types  of 
approach  are  clearly  complementary  at  this  point.  However,  the  assumption  of 
a  slowly  varying  wave  train  is  unable  to  cope  with  a  situation  in  which  two 
quite  different  Fourier  components  of  the  surface  wave  are  involved.  In  the 


40.  L.M.  Brekhovskikh,  et  al,  "On  the  resonant  generation  of  internal  waves 
through  nonlinear  interactions  of  surface  waves,"  Bull.  (Dzv.)  Acad.  Sci. 
USSR,  Atmos,  and  Oceanic  Phys.,  8,  1972,  192. 
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same  way,  the  simple  resonance  theory  is  not  well  suited  to  consider  the 


blockage  phenomenon  that  occurs  when  C  -  c  is  not  small  and  the  adverse 

g 

current  is  substantial. 


3.  Ko's  Approach 

The  Ko  approach  to  the  interaction  of  internal  waves  and  surface  gravity 
waves  is  given  in  Ref  24.  He  defines  the  internal  wave  induced  fractional 
changes  in  surface  wave  amplitude  a  and  slope  m  as  follows: 


d  *  3 

max  min 

a  +  a  . 

max  min 


and 


m* 


m  -  m  . 
max  min 

m  +  m  . 
max  min 


(97) 


(98) 


Instead  of  a  general  solution  to  Eq  (1)  and  (9)  for  this  problem,  this  study 
was  directed  toward  the  one-dimensional  wave  interaction  corresponding  to 
their  tank  experiment.  The  surface  current  U  from  the  internal  wave  is- 
generated  by  the  interfacial  waves  of  a  two-fluid  system  with  frequency  tl  and 
phase  speed  C  for  time  t  >  0.  An  independent  surface  wave  generator  with 
frequency  u)  is  used  to  produce  for  all  time  a  single  surface  wave  train  with 
wavenumber  k^  and  unperturbed  group  speed  The  problem  was  to  determine 

the  change  in  the  surface  wave  characteristics  as  functions  of  both  position  x 
along  the  tank  and  time  t,  when  the  internal  wave  induced  surface  current 
field  U  was  present. 

For  these  one-dimensional  waves,  the  governing  relations,  ie,  Eq  (1)  and 
(7),  were  simplified  as  follows.  Equation  (1)  became 


at  *  *  ^^ax  "  "Sx  ’ 

where  =  %(g/k)^  and  U  =  U(x,t).  Equation  (7)  reduced  to 


^  +  h  [(C.  +  U)A2]  =  -^2 


at  ax  g 


(100) 


which  was  written  in  the  alternative  form 


Mi  +  (e  +  U)^  =  -^a2^  -  a2^ 
at  g  ^■'ax  2  ax  ax  ' 


(101) 


Both  Eq  (99)  and  (101)  were  written  in  a  characteristic  form  for  which  the 
solution  could  be  readily  obtained.  These  characteristics  were  defined  by 


TT  =  c  +  U  , 
dt  g  ’ 


(102) 


for  which  the  two  governing  equations  reduce  to 
dt  ~  ‘^ax 


(103) 


-  .  1a2  §0  .  a2 _ & 

2  ax  ax  • 


(104) 


In  other  words,  in  the  characteristic  coordinate  system  defined  by  integrating 
Eq  (102)  such  that  x  =  x(t,4)*  t  =  t,  where  4  =  constant  denotes  a  charac¬ 
teristic,  the  governing  equations  [Eq  (103)  and  (104)]  become 


aAnk  _  au 

at  ■  "ax 

i  t 


(105) 


aEnA^  _  3  ay 
at  ■  ■  2  ax 

i  t 


(106) 


To  simulate  the  tank  experiment  the  surface  current  field  was  represented  by 


(107) 


U  =  sin  K(x  -  Ct)H(x)H(t  -  x/C)  , 
where  H(x)  represents  the  Heaviside  function. 

To  simplify  the  analysis  further,  it  was  assumed  that  the  magnitude  of 
the  surface  current  u^  was  small  compared  with  the  minimum  of  the  group  speed 
for  surface  waves,  ie. 


where  c^^  =  %(g/k^)^,  the  unperturbed  group  speed.  Then  it  was  assumed  that 
the  solutions  could  be  represented  by  the  following  perturbation  series  in  e: 


X  =  x^  +  exj  +  e^X2  +  . . .  , 

£nk  =  Ank^  +  e£nki  +  e^£nk2  +  . . .  , 

inA^  =  AnA^  +  eitnAf  +  e^£nA|  + _  , 


where  c  =  -^c  *£nki. 
gi  go  ^ 

Substituting  Eq  (109)  into  Eq  (102),  (105),  and  (106)  and  grouping  terms 
of  the  same  order  in  e  gave  for  0(e®): 


3x  dk 

_ o  _  _ o 

8t  ■  ^go  ’  at 


aA2 
_ o 

at 


ac 


=  -A 


2_8o  = 

o  ax 


0  . 


(110) 


These  equations  describe  the  unperturbed  surface  wave  field,  which  was  assumed 
to  be  a  uniformly  generated  monochromatic  wave  train.  Integration  of  Eq  (110) 
gives 


k  =  const.,  a2  =  const.,  x  =  c  (x  ~  4)  , 
o  ’  o  o  go  ^  ’ 


(111) 
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where  4  denotes  the  intersection  of  the  characteristic  with  the  t  axis.  For 
0(e)  the  governing  equations  become 


=  u  [K  cos  K(x  -  Ct)H(x  )H(t  -  x  /C) 

Ol  o  o  o  o 

+  sin  K(x^  -  Ct)6(x^)H(t  -  x^/C) 

-  sin  K(x  -  Ct)H(x  )6(t  -  x  /C)l  ,  (112) 

o  o  o 


3£nA?  _  3  3Anki 
3t  “2  3x 


(113) 


and 

e-^  =  u  sin  K(x  -  Ct)H(x  )H(t  -  x  /C)  +  ec  ,  (114) 

Ol  o  o  o  o  Si 

where  eAnk^  and  e£nAf  represent  the  first-order  change  of  the  wavenumber  and 

the  square  of  the  amplitude,  respectively,  while  ex^  gives  the  deviation  from 
the  unperturbed  characteristic. 

Using  the  zeroth-order  characteristic  gives 


X  -  Ct 
o 


c)t 


(115) 


Then  the  equations  could  be  integrated  from  t  =  0  to  x  along  a  given 
characteristic  4  =  const.  The  solution  to  Eq  (112)  was  obtained  first,  then 
used  on  the  right-hand  sides  of  Eq  (113)  and  (114).  The  solution  for  the 
wavenumber  is  given  by 


eAnkj  = 

o 


u 

- 7;H(x  )  sin  K(x  -  Ct)H(t  -  x  /C) 

C  •  L  O  O  O 

go  L 


and  for  the  amplitude  by 


-^Kx  cos  KC  cos  KC 
c  o 

go 


&  •  ■)]  ■(■  -  y 


(117) 


For  k  =  +  Ak  and  A  =  +  M,  the  right-hand  sides  of  Eq  (116)  and  (117) 

give  the  changes  in  wavenvunber  (Ak/k^)  and  wave  amplitude  (AA/A^), 
respectively,  for  small  changes.  In  practice,  their  interest  was  in  the 
change  of  the  local  surface  wave  slope,  which  is  related  to  the  changes  in 
wavenumber  and  amplitude  by 


/>  ™  /I  h  «  A 
An —  =  £n; —  +  En^—  , 
m  k  A  ’ 

o  o  o 

or  for  small  variations  by 
Am  _  AA  ^  Ak 


(118) 


(119) 


They  found  it  interesting  to  also  calculate  the  maximum  changes  in  the 
surface  wavenumber,  the  amplitude  and  the  slope  at  a  given  x-station,  as  well 
as  the  corresponding  phase  of  the  interval  wave  current  where  the  maximum 
effect  occurs.  For  the  region  where  both  Heaviside  functions  in  Eq  (116)  and 
(117)  are  equal  to  1,  it  is  appropriate  to  drop  the  Heaviside  function 
notation.  Furthermore,  they  define 

(|i  =  K(x^  -  Ct)  (120) 

and 
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Then  Eq  (116),  (117)  and  (118)  can  be  rewritten  as 


(123 


(124 


(125 


For  a  given  and  C/c  ,  the  maxiinum  or  minimum  change  occurs  at  the 
o 


phase  of  the  internal  wave  where  8()/9(b  =  0;  the  quantity  inside  the 
parentheses  can  be  any  one  of  the  three  quantities  given  on  the  left  of  Eq 
(123),  (124),  and  (125).  The  results  of  performing  the  indicated 


differentiations  are 


-  C  c 


cos  B  +  - ^Kx  sin  B 

c  o 
go 


sin  B  -  - ^Kx  cos  B 

c  o 
go 


C  c 


— cos  B  +  -^1 


Kx  sin  B 


:: —  I  I  °  - n  I  sin  B  -  - ^Kx  cos  B 

cll  c-CI  c  o 

go  L\  go  /  go 


The  corresponding  maximuoi  or  minimum  values  are  obtained  by  substituting  Eq 
(126),  (127),  and  (128)  into  Eq  (123),  (124),  and  (125),  respectively. 

Since  the  solutions  of  Eq  (116)  and  (117)  appear  to  be  inversely 
proportional  to  C  -  c^^,  a  singular  behavior  may  seem  probable  at  resonance. 
However,  with  a  proper  limiting  procedure,  it  can  be  shown  that,  for  c^^  C 


£n^  = 
k 

o 


jr-H(x  )H(t  -  X  /C)[sin  K(x  -  Ct)  +  Kx  cos  K(x  -  Ct)]  (129) 

O  O  O  O  O 


£nf-  =  -  77H(x  )H(t  -  X  /C) I (4  -  K^x^)  sin  K(x  -  Ct) 

n  O  O  O  O 

O 


+  5Kx  cos  K(x  -  Ct)]  . 
o  o 


(130) 


Notice  the  strong  dependence  of  the  solution  on  kx^,  which  effectively 
measures  the  distance  or  "time”  of  interaction.  Further,  in  this  model,  the 
energy  transfers  only  from  the  internal  wave,  and  there  is  no  provision  to 
account  for  the  limited  energy  available  in  the  internal  wave.  Without 


considering  the  other  nonlinear  effects  and  the  energy  loss  of  the  internal 


wave  train,  the  linearized  solutions  for  the  resonance  case  predict  a 
monotonic  increase  of  the  effects  with  distance  as  shown  by  Eq  (129)  and 

(130). 


Both  the  theoretical  and  experimental  results  for  slope  and  amplitude 
modulations  show  that,  when  normalized  by  the  ratio  of  induced  surface  current 
to  internal  wave  phase  speed  (ie,  u^/C),  the  "resonant"  case  (ie,  C  -  c^^) 
modulation  magnitudes  are  dependent  only  on  the  interaction  "distance"  kx  and 
are  independent  of  the  magnitudes  of  the  other  physical  parameters.  The 
maximum  normalized  (by  u^/C)  slope  and  amplitude  interaction  effects  are  found 
to  occur,  for  both  theory  and  experiment,  when  the  phase  speed  of  the  internal 
wave  and  the  group  speed  of  the  surface  wave  are  matched  (ie,  C  =  c  ). 


4.  HUGHES'  APPROACH 

Hughes  (Ref  17)  begins  with  Eq  (1),  (2),  and  (5),  which  in  his  notation 
reduce  to 


*  o)  =  0 


ko  =  const. 


8A  8  , .  ,  .  8u» 

^  *  ajIACu  ♦  'g)l  =  A.  ,  Cg  =  5J- 


Note  that  he  has  added  a  source  (or  "growth/decay")  term,  which, 
n 


(131) 

(132) 

(133) 


we  shall  see, 


leads  to  Z  S.  of  Eq  (19).  By  manipulation  of  his  variable  dependence,  he 
i=l  ^ 

arrives  at  Eq  (8)  in  the  form 


where  he  uses 


3 

A  =  1  S  =  aCki.ka)  +  p(ki,k2)N  -  y(ki,k2)N2  (135) 

i=l 

in  place  of  the  relations  in  Section  II. £.  He  justifies  ignoring  a(ki,k2)  and 
observes  that  the  third  term  represents  a  nonlinear  limitation  to  the 
exponential  growth,  and  while  it  does  not  allow  for  overshoot  effects,  it  does 
provide  an  inherent  upper  limit  to  the  spectral  growth. 

Let  k  be  the  initial  value  of  ki  in  the  sense  that  ki(k  ,x,t)  =  k  at  x 

O  *  ^0^0 

=  t  =  0.  If  the  wave  spectrum  in  the  absence  of  interaction  is  stationary  or 
growing/decaying  very  slowly,  the  left  side  of  Eq  (134)  is  zero;  therefore  at 
equilibrium,  A  =  0,  so  that 


Y(k^,k2)  =  p(k^,k2)/N^  , 


(136) 


where  N^(k^,k2)  is  the  unperturbed  equilibrium  spectral  value.  It  is  also 
assumed  that  this  condition  holds  for  perturbed  wave  numbers,  ie,  during 
surface-internal  wave  interaction,  but  with  k^  replaced  by  ki.  This  is 
equivalent  to  assuming  that  the  wind  shifts  the  spectral  amplitude  for  given  k 
toward  the  equilibrium  spectral  amplitude  of  the  same  wavenumber.  Under  these 
conditions,  Eq  (134)  becomes 


|s  ♦  (U(.,t)  »  c^l|S  =  P(k,,k2)»  -  jTjl^ 


(137) 


To  complete  the  mathematical  statement  of  the  problem,  ie,  Eq  (131),  (132), 
and  (136),  initial  conditions  on  k^  and  N  are  imposed: 


ki(t  =  0,x,y)  =  k^ 

N(t  =  0,x,y)  =  N^(k^,k2)  .  (138) 
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Using  the  method  of  characteristics  on  Eq  (131),  can  be  obtained  by 
solving  the  pair  of  coupled  nonlinear  first-order  ordinary  differential 
equations 

^  =  U(x.t)  .  c^ 


dt  '^^ax 


(139) 


x=x(t) 


where  x(0)  is  arbitrary.  Along  the  path  x(t)  with  ki(t)  known  (in  principle), 
Eq  (137)  becomes 


dt 


=  pIki(t),k2lN  -  p2N2/N^(kj(t),k2)  , 


(140) 


which  is  just  a  particular  form  of  Eq  (9).  Equation  (140)  is  a  simplified 
form  of  Eq  (19),  whose  solution,  after  integrating  once  by  parts  on  r)  and 
using  Eq  (139),  is 


N(k,(t),k2)  =  Njki(t),k2l  j  1  -  [ki(t),k2l  J  ^ 

o 


dt|i 

ak 


:(«)  “'(-j 

^  'x(n)  \  n 


exp  I  -  j  P(ki(q),k2)dq  I  dn 

n 


-1 


(141) 


If  the  orbital  velocity  of  the  internal  waves  is  represented  by  U  =  U(x 
Ct),  where  C  is  the  phase  velocity  of  the  internal  wave,  then  trapped  waves 
exist  between  locations  where 


c  +  U(x  -  Ct)  -  C  =  0 
8 


(142) 


as  described  in  Ref  38  and  in  Ref  1,  p  78-81.  Hughes  (Ref  17)  discussed 
some  properties  of  the  solution  of  Eq  (141),  including  the  trapped  waves. 
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Hughes  (Ref  17)  aext  considers  only  steady,  dispersionless  current 
fields.  Therefore,  he  particularizes  the  theory  to  the  case  of  U  =  U(x  -  Ct) 
and  |U|/C  -►  0.  If  all  powers  of  |Ul/C  except  the  first  are  ignored,  Eq  (139) 
and  (141)  have  the  solution 


X  =  X  +  c  t  +  0(U) 
o  go 


(143a) 


-  /  U'[x^  +  (c^^  -  Onldn^ 


(143b) 


dN 


N(k,,k2)  =  N^(ki,k2)  +  ki^  J  U’[x^  -  (Cg  -  On] 

O 

•  exp  [-P(ki,k2)(t  -  n)]dn. 

■  (A  ■ 


(143c) 


where  U'  =  and  c 


In  terms  of  the  wave  displacement  spectrum  ijj,  with  U  expressed  as  a  sum  of 
different  internal  wave  wavenumber  components,  ie. 


N 


U(x 


,t)  =  B  sin  [K  (x  -  Ct)  +  0  ]  , 

Z  j  n  n  n 


(144) 


n=l 


the  integration  can  be  carried  out  to  yield 


"  B  ak,  an 

«l»(ki,k2)  =  4<o(ki,k2)  +  Yj  "V”  ak7 

n=l 


0/CK^)C^(x,t)  +  (c^/C  -  l)D^(x,t) 
(P/cK  )2  +  (c„/C  -  1)^ 


(145) 


where  for  t  >  30  min 


C  (7<,t)  =  cos  [K  (x  -  Ct)  +  6  ] 
n  n  n 


D  (x,t)  =  sin  (K  (x  -  Ct)  +  0  ]. 
n  ’  n  n 


(146) 


In  Eq  (145)  note  that  exact  resonance  occurs  when  c  =  C  and  each  component  of 

S 

the  surface  wave  field  is  additively  perturbed  by  each  component  of  the 
internal  wave  field  (in  this  approximation). 

From  ordinary  resonance  analysis,  the  phase  angle  for  each 
perturbation  component  is  given  by 


%  =  c  /C  -  i 
8 


(147) 


and  thus  for  weak  damping,  increases  from  -180°  to  0°  as  c^  increases  from 
much  less  than  C  to  much  greater  than  C.  As  ^/Ck^  0,  the  range  of  c^  over 

which  changes  appreciably  becomes  more  closely  confined  to  the  resonance 
point  c^  =  C.  The  phenomenon  of  trapping  is  degenerate  in  this  linear 
approximation,  and  only  those  waves  that  satisfy  the  resonance  condition  [Eq 
(142)  with  U/C  ■*  0]  can  be  said  to  be  "trapped." 

With  this  description  and  with  a  least-squares  regression,  the 
sensitivity  and  phase  of  the  mean  square  slope  can  be  defined.  Let  s^  repre¬ 
sent  the  total  surface  slope  variance  (the  sura  of  the  variances  of  the  two 
slope  components).  Then 


^  =  SS  (1^1  +  k|)»l>(ki,k2)dkidk2 
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and  from  Eq  (145) 


f4 


^  -  T^\  1 


- —  1=  y  B  R  sin  (k  (x  -  Ct)  +  0  +  |J  1 

~2“  /  ^  n  n  ‘  n  n  '^n 

/  n=l 


(149) 


N 


=  ^  B  R  islnfk  (x  -  Ct)  +  6  ]  cos  M 

/  ■  n  n  (  n  n  n 


n=l 


+  cos  [k  (x  -  Ct)  +  6  ]  sin  |J  > 
n  n  n  I 


where  R^  and  are,  respectively,  the  sensitivity  and  phase  of  the  component 


s^  with  wavenumber  k  ,  ie. 
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In  accordance  with  the  measurement  process,  an  overall  sensitivity  R  and 
phase  p  may  be  determined  by  correlating  the  perturbations  in  ^  with  phase- 
shifted  current  values  and  defining  p  as  the  phase  shift  that  maximizes  the 
correlation  and  R  as  the  corresponding  regression  line  slope.  By  performing 
the  required  manipulations  to  Eq  (144)  and  (139),  it  can  be  shown  that 
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Values  for  and  |j^  are  obtained  by  evaluating  the  integrals  given  in  Eq 
(150),  but  before  this  can  proceed  it  is  necessary  to  estimate  4<Q(l^i>lt2) 

I 

each  data  set. 

Measured  values  of  4*  integrated  in  angle  at  various  values  of  x,t  within 

the  internal  wave  field  can  also  be  obtained.  That  is,  with  ki  =  k  cos  0,  kg 
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=  k  sin  0  and  H(k)  =  /  d0tl)(k,0). 
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where 


P  =  a[(u^cos  0/c]{O.Ol  +  0.016*|cos  0|*Ujj./c} 


•{1  -  exp[-8.9(u^,/c  -  0.03)^]} 


(153) 


and  0  is  the  angle  between  the  surface  wave  propagation  and  the  wind 
direction. 

The  theoretical  and  measured  values  of  the  sensitivity  gave  good 
agreement,  while  for  the  phase  p  the  agreement  was  poor  except  that  all 
theoretical  and  measured  values  lay  in  the  same  quadrant.  An  illustration  of 
the  agreement  between  the  theoretical  and  measured  forms  for  the  integrated 
spectrum  H(k)  is  given  in  Fig  4  of  Ref  17.  The  agreement  was  satisfactory 
except  at  high  wavenumbers. 
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F.  WAVE  ATTENUATION 


Wave  attenuation  coefficients  are  found  in  Ref  1,  p  51-54.  The  simplest 
case  is  that  of  deep  water  and  a  clean  surface.  Although  the  viscous  stresses 
in  the  surface  layer  are  important,  the  layer  itself  makes  a  negligible 
contribution  to  the  overall  rate  of  energy  dissipation  since  the  rates  of 
strain  there  are  no  larger  than  in  the  irrotational  flow  and  the  layer 
thickness  is  small.  The  energy  losses  arise  almost  entirely  from  the 
straining  of  the  irrotational  motion,  so  that  the  attenuation  coefficient  is 

Yy  =  2vk2  ,  (154) 

where  the  kinematic  viscosity  of  water  is  V  -  0.01  cm  ^sec  ^  at  20“C  and 
35‘*/oo  salinity.  The  energy  density  of  the  wave  field  decreases  as  exp(-2Y^t) 
and  the  amplitude  as  expC-y^t). 

With  a  densely  packed  surface  film,  the  attenuation  factor  Yf  resulting 
from  the  surface  layer  alone  is  now 

Yf  =  ^vkVa/2v  (155) 

in  deep  water.  The  variation  of  y^  and  Y^  with  wavelength  and  frequency  is 
illustrated  in  Fig  3.1  of  Ref  1.  It  is  evident  that,  except  for  the  very 
short  capillaries,  Y^  Yy-  Short  waves  incident  upon  a  slick  are  quite 
rapidly  attenuated  and,  since  the  apparent  smoothness  of  the  sea  surface  is 
very  much  dependent  on  small-scale  components,  this  attenuation  evidently 
accounts  for  the  smooth  appearance  of  oil  slicks  in  the  ocean. 
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IV.  PRELIMINARY  CONCLUSIONS  AND  RECOMMENDATIONS 


The  purpose  of  this  study  of  the  state  of  the  art  of  sea  surface  modeling 
was  to  determine  the  feasibility  of  designing  a  software  model  that  could 
adequately  describe  the  sea  surface  for  treating  (1)  acoustic  torpedo  surface 
reverberation  in  shallow  water  and  (2)  remote  radar  surface  imaging.  From  the 
above  material  we  come  to  the  preliminary  conclusion  that  such  a  sea  surface 
software  model  is  presently  possible.  We  now  outline  how  such  a  software 
model  might  be  developed  and  implemented. 

The  basic  tool  for  the  sea  surface  model  would  be  the  radiative  transfer 
equation,  Eq  (19),  as  described  in  Section  II. F.  Since,  as  seen  in  Section 
III,  the  interaction  of  surface  and  internal  waves  is  the  best-understood 
phenomenon,  this  provides  a  good  area  in  which  to  start  developing  our  sea 
surface  software  model.  First,  the  radiative  transfer  relation  of  Ref  17 
would  be  implemented  based  on  the  information  contained  in  Section  III.E.4. 
This  would  be  compared  with  the  experimental  results  of  Ref  39  (Section 
III.E.4)  and  Ref  24  (Section  III.E.3).  Then  additional  radiative  transfer 
source  terms  from  Section  II.F  would  be  included  and  the  experimental 
comparison  would  be  repeated.  As  a  further  check,  the  model  would  be 
reconciled  with  Eq  (88)  and  (93) . 

Next  the  model  would  be  extended  to  agree  with  the  JONSWAP  spectrum  of  Eq 
(27)  and  the  directional  relation  of  Eq  (17).  Then  the  model  would  be 
generalized  to  consider  variable  surface  currents  and  checked  against  Eq  (41), 
(43),  (49),  and  (58).  Finally  swell  conditions  would  be  added  and  tests  made 
against  Eq  (73).  At  this  point  all  four  phenomena  of  interest  (ie,  wind 
velocity,  surface  currents,  swell,  and  internal  waves)  would  be  present  in  the 
development . 


To  be  useful  for  the  stated  naval  system  problem  modeling,  the  sea 
surface  software  model  should  be  treated  to  yield  the  spectra  for  slopes, 
velocities,  and  accelerations  and  then  extensively  tested  by  varying  the  many 
variables  and  checking  the  results.  The  material  following  Eq  (148)  presents 
an  interesting  area  of  experimental  verification  for  slope  properties,  as  does 
the  information  in  Ref  24,  Section  III.E.3. 

The  total  software  package  for  the  statistical  modeling  of  the  sea 
surface  would  then  be  ready  for  application  to  modeling  the  acoustic  torpedo 
surface  reverberation  in  shallow  water  and  remote  radar  sea  surface  imaging. 
This  would  be  accomplished  by  formulating  (1)  software  that  gives  a 
statistical  description  of  surface  reverberation  and  (2)  software  that  gives  a 
statistical  description  of  imaging  for  remote  radar  sensing.  The  latter 
should  be  extended  to  address  the  phenomena  of  the  directional  attenuation  of 
surface  wave  trains  and  solitary  waves. 


REFERENCES 


1.  O.M.  Phillips,  The  Dynamics  of  the  Upper  Ocean,  Cambridge  University 
Press,  1977. 

2.  F.  Unsell,  "Steady  wave  patterns  on  a  non-uniform  fluid  flow," 

J.  Fluid  Mech..  9,  1960,  333-346. 

3.  G.B.  Whitham,  "A  note  on  group  velocity,"  J.  Fluid  Mech. .  9,  1960, 
347-352. 

4.  F.P.  Bretherton  and  C.J.R.  Garrett,  "Wavetrains  in  inhomogeneous  moving 
media,"  Proc.  Roy.  Soc.,  A,  302,  1969,  529. 

5.  F.P.  Bretherton,  "The  general  linearized  theory  of  wave  propagation,"  in 
Mathematical  Problems  in  the  Geophysical  Sciences,  vol.  1,  61-102,  Amer. 

Nath.  Soc.,  Providence,  RI ,  1971. 

6.  M.S.  Longuet-Higgins  and  R.W.  Stewart,  "Changes  in  the  form  of  short 
gravity  waves  on  long  waves  in  tidal  currents,"  J.  Fluid  Mech. ,  8,  1960, 
565-583. 

7.  M.S.  Longuet-Higgins  and  R.W.  Stewart,  "The  changes  in  amplitude  of  short 
gravity  waves  on  steady  non-uniform  currents,"  J.  Fluid  Mech.,  10,  1961, 


529-549. 


8.  6.B.  Whitham,  "Mass  momentuni  and  energy  flux  in  water  waves," 

J.  Fluid  Mech. .  12,  1962,  135-147. 

9.  T.P.  Barnett  and  K.E.  Kenyon,  "Recent  advances  in  the  study  of  wind 
waves,"  Rep .  Prog .  Phys . ,  38,  1975,  667-729. 

10.  H.  Mitsuyasu,  "Measurement  of  the  high-frequency  spectrum  of  ocean 
surface  waves,"  J .  Phys .  Ocean . .  7,  1977,  882-891. 

11.  S.A.  Kitaigorodskii ,  "The  statistical  characteristics  of  wind-generated 
short  gravity  waves,"  in  Spaceborne  Synthetic  Aperture  Radar  for  Oceanograph 
R.C.  Beal,  et  al,  eds..  The  Johns  Hopkins  University  Press,  1981. 

12.  C.S.  Cote,  et  al.  Meteorological  Papers,  2(6),  W.J.  Pierson,  Jr.,  ed., 
1960. 


13.  M.S.  Longuet-Higgins ,  D.E.  Cartwright,  and  N.D.  Smith,  Ocean  Wave 
Spectra.  Prentice-Hall,  1963,  111-136. 

14-  D.E.  Cartwright  and  N.D.  Smith,  Buoy  Technology,  Marine  Technology 
Society,  Washington,  D.C.,  1964,  112-121. 

15.  J.A,  Ewing,  J.  Marine  Res.,  27,  1969,  163-171. 

16.  O.M.  Phillips,  "The  structure  of  short  gravity  waves  on  the  ocean 
surface,"  in  Spaceborne  Synthetic  Aperture  Radar  for  Oceanography,  R.C.  Beal 


et  al,  eds..  The  Johns  Hopkins  University  Press,  1981. 


17.  B.A.  Hughes,  The  effect  of  internal  waves  on  surface  wind  waves.  2: 
Theoretical  analysis,"  J.  Geophys.  Res.,  83 ,  1978,  455-465. 

18.  H.  Hasselmann,  Basic  Developments  in  Fluid  Dynamics,  Vol .  2,  M.  Holt, 
ed..  Academic  Press,  1968,  117-182. 

19.  M.S.  Longuet-Higgins ,  "On  the  nonlinear  transfer  of  energy  in  the  peak  of 
a  gravity-wave  spectrum:  a  simplified  model,"  Proc.  Roy.  Soc.,  A,  437 ,  1976, 
311-328. 

20.  T.P.  Barnett,  "On  the  generation,  description,  and  prediction  of  ocean 
wind  waves,"  Ph.D.  dissertation.  University  of  California,  San  Diego,  1966. 

21.  T.P.  Barnett,  "On  the  generation,  description,  and  prediction  of  ocean 
wind  waves,"  J .  Geophys .  Res . ,  73,  1968,  513-529. 

22.  J.A.  Ewing,  "A  numerical  wave  prediction  for  the  North  Atlantic  Ocean," 
Deutsche  Hydrog.  Z. ,  24 ,  1971,  241-261. 

23.  O.M.  Phillips,  "On  the  interactions  between  internal  and  surface  waves," 
Phys .  Atmos .  and  Oceans ,  9,  1973,  954-961. 

24.  J.E.  Lewis,  B.M.  Lake,  and  D.R.S.  Ko,  "On  the  interaction  of  internal 
waves  and  surface  gravity  waves,"  J.  Fluid  Mech. ,  ^(4),  1974,  773-800. 


25.  C.  Garrett  and  J.  Smith,  "On  the  interaction  between  long  and  short 
surface  waves,"  J.  Phys.  Oceanogr.,  6,  1976,  925-930. 


26.  O.H.  Phillips,  "The  dispersion  of  short  wavelets  in  the  presence  of  a 
dominant  long  wave,"  J.  Fluid  Mech.,  107,  1981,  465-485. 

27.  O.H.  Shemdin,  "Modulation  of  centimetric  waves  by  long  gravity  waves: 
Progress  report  on  field  and  laboratory  results,"  NATO  Conference  on  Turbulent 
Fluxes  Through  the  Sea  Surface,  Wave  Dynamics  and  Prediction,  A.  Farve  and  K. 
Hasselmann,  eds.,  NATO  Conference  Series  V,  vol.  1. 

28.  E.D.  Graves  and  O.H.  Shemdin,  "An  investigation  of  the  modulation  of 
capillary  and  short  gravity  waves  in  the  open  ocean,"  J.  Geophys.  Res.. 

85(C9),  1980,  5019-5024. 

29.  G.R.  Valenzuela  and  J.W.  Wright,  "Modulation  of  short  gravity-capillary 
waves  by  longer-scale  periodic  flows — a  higher-order  theory,"  Radio  Science, 
14(6),  1979,  1099-1110. 

30.  M.J.H.  Fox,  "On  the  nonlinear  transfer  of  energy  in  the  peak  of  a 
gravity-wave  spectrum — I,"  Proc.  Roy.  Soc..  A  348,  1976,  467-483. 

31.  K.  Hasselmann,  et  al,  "Measurements  of  wind  wave  growth  and  swell  decay 
during  the  Joint  North  Sea  Wave  Project  (JONSWAP),"  Herausgegeben  von  Deutsche 
Hydrograph.  Institut,  Reihe  A,  11,  1973. 

32.  D.  Holliday,  "Nonlinear  gravity-capillary  surface  waves  in  a  slowly 
varying  current,"  J.  Fluid  Mech.,  5^(4),  1973,  797-802. 


33.  M.S.  Longuet-Higgins ,  "The  instabilities  of  gravity  waves  of  finite 
amplitude  in  deep  water,  Part  I:  Superharmonics,"  Proc.  Roy.  Soc.  A  360, 

1978,  471-488. 

34.  R.C.  Beal,  "Spaceborne  imaging  radar:  monitoring  of  ocean  waves," 
Science.  208,  1980,  1373-1375. 

35.  D.D.  Evans  and  D.H.  Shemdin,  "An  investigation  of  the  modulation  of 
capillary  and  short  gravity  waves  in  the  open  ocean,"  J .  Geophys .  Res . .  85, 
1980,  5019-5024. 

36.  H.  Mitsuyasu,  "Interactions  between  water  waves  and  wind  (1)," 

Rep.  Res.  Appl.  Mech.,  Kyushu  University,  1^,  1966,  67-88. 

37.  O.M.  Phillips  and  M.L.  Banner,  "Wave  breaking  in  the  presence  of  wind 

drift  and  swell,"  J.  Fluid  Mech. .  1974,  625-640. 

38.  A.E.  Gargett  and  B.A.  Hughes,  "On  the  interaction  of  surface  and  internal 
waves,"  J.  Fluid  Mech. .  52(1),  1972,  179-191. 

39.  B.A.  Hughes  and  H.L.  Grant,  "The  effect  of  internal  waves  on  surface  wind 
waves.  1:  Experimental  measurements,"  J.  Geophys.  Res.,  83 ,  1978,  443-454. 

40.  L.M.  Brekhovskikh,  et  al,  "On  the  resonant  generation  of  internal  waves 
through  nonlinear  interactions  of  surface  waves,"  Bull.  (Dzv.)  Acad.  Sci. 

USSR,  Atmos,  and  Oceanic  Phys.,  8,  1972,  192. 


Oceanography,  R.C.  Beal,  et  al,  eds.,  Johns  Hopkins  University  Press,  1981. 


42.  R.O.  Hagar,  "SAR  ocean  imaging  mechanics,"  in  Spaceborne  Synthetic 
Aperture  Radar  for  Oceanography,  R.C.  Beal,  et  al,  eds.,  Johns  Hopkins 
University  Press,  1981. 

43.  G.R.  Valenzuela,  "An  asymptotic  formulation  of  SAR  images  of  the 
dynamical  ocean  surfaces,"  Radio  Science.  15(1),  1980,  105-115. 


44.  W.R.  Alpers,  O.B.  Ross,  and  C.L.  Rufenach,  "On  the  detectability  of  ocean 
surface  waves  by  real  and  synthetic  aperture  radar,"  J.  Geophys.  Res.,  86(C7) , 
1981,  6481-6498. 


APPENDIX:  IMAGING  FOR  REMOTE  RADAR  SENSING 


Beal  (Ref  41)  discusses  the  fundamentals  of  aperture  synthesis  for  the 
Seasat  SAR  (synthetic  aperture  radar).  More  detailed  discussions  of  imaging 
for  remote  radar  sensing  are  found  in  Ref  42-44.  Seasat  SAR  obtains  its  high 
resolution  by  artificially  synthesizing  an  aperture  many  kilometers  long  in 
space.  For  studying  the  ocean,  SAR  has  two  essential  advantages  over  other 
high-resolution  remote  sensors:  it  collects  its  data  through  cloud  cover  and 
without  help  from  the  sun.  It  has  demonstrated  the  potential  ability  to  image 
and  globally  monitor  ocean  wave  systems.  SAR  is  sensitive  primarily  to  the 
structure  of  the  short  gravity  waves  (30-cm  waves  for  Seasat  SAR)  on  the  ocean 
surface  that  were  shown  in  Section  III  to  be  sensitive  to  wind  velocity, 
surface  currents,  swell,  and  internal  waves. 

According  to  classical  diffraction  theory,  the  angular  resolution  of  any 
transmitting  or  receiving  system,  including  SAR,  is  ultimately  limited  by  the 
size  of  its  aperture.  Expressed  simply 

♦  =  1/2  ,  (A-1) 
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where  <ti  is  the  angular  resolution  and  £  is  the  size  of  the  aperture  expressed 
in  wavelengths.  When  the  same  aperture  is  used  for  both  the  transmitter  and 
the  receiver  of  energy,  as  it  is  in  SAR,  the  resulting  angular  resolution  ({>' 
is  effectively  halved  (ie,  <|>'  =  (}»/2).  The  Seasat  orbited  at  a  slant  range  h 
(20°  from  its  nadir)  of  about  10^  m  and  contained  a  SAR  that  operated  at  a 
radar  wavelength  of  0.2  m.  The  length  L  of  the  synthetic  aperture  is  related 
to  the  azimuth  resolution  r  (the  ground  resolution  in  the  flight  direction  of 
the  satellite)  by 

r  =  (|('h  =  ^h  .  (A-2) 

d  la 

An  azimuth  resolution  of  25  m  requires  a  synthetic  aperture  of  4  km.  A 
satellite  traveling  at  an  altitude  of  about  10^  m  has  an  orbital  velocity  of 
about  8  km/s  and,  therefore,  requires  about  0.5  s  to  synthesize  the  required 
aperture . 

Note  that  the  length  of  the  real  aperture  D  on  the  spacecraft  has  not 
explicitly  entered  into  the  equations.  However,  the  real  aperture  must  be 
short  enough  to  allow  a  particular  point  on  the  ground  to  remain  entirely 
within  the  real  beam  0  =  \/D  during  the  aperture  synthesis  interval.  This 
leads  to  the  relationship 

L  =  hg  .  (A-3 

Combining  Eq  (A-2)  and  (A-3)  leads  to  the  fundamental  lower  limit  for 


resolution 


Therefore,  if  a  ground  resolution  of  25  m  is  desired,  the  real  aperture  can  be 
no  longer  than  50  m  to  maintain  a  point  in  the  real  beam  for  a  sufficiently 
long  time.  The  resolution  limit  is  independent  of  range  because  the  time 
during  which  a  particular  point  is  illuminated  increases  with  range,  thus 
allowing  a  corresponding  larger  synthetic  aperture  to  be  formed.  The  Seasat 
SAR  has  a  real  aperture  length  of  about  12  m  and  was,  therefore,  theoretically 
capable  of  an  azimuth  resolution  of  6  m.  This  assumes  that  the  maximum 
allowable  synthetic  aperture  of  16  km  could  be  constructed,  a  distance  that  is 
covered  by  the  satellite  in  about  2  s. 

Although  a  resolution  of  6  m  in  azimuth  is  theoretically  possible  from 
Seasat,  the  imagery  is  typically  processed  to  yield  only  25  m.  That  is,  only 
4  km  of  synthetic  aperture  (about  0.5  s  of  data)  are  simultaneously  processed. 
Any  4  m  from  the  total  16  km  length  will  satisfactorily  produce  a  resolution 
of  25  m.  Moreover,  since  the  predominant  system  noise  in  the  SAR  is  caused  by 
"coherent  speckle,"  which  tends  to  be  Rayleigh  distributed  in  amplitude  (ie, 
there  are  wide  variations  in  the  reflected  signal  from  a  resolution  element  as 
the  illumination  angle  is  slightly  varied),  each  4-km  segment  potentially 
allows  an  independent  sample  of  the  "average"  reflectivity  distribution. 
Therefore,  the  variance  within  a  scene  can  be  reduced  considerably  by 
separately  processing  and  combining  four  independent  25-m  images. 

Thus  far  the  discussion  has  dealt  only  with  the  very  basic  criteria  for 
obtaining  high  azimuth  resolution  along  the  velocity  vector  of  the  satellite - 
We  have  not  examined  the  actual  mechanics  of  collecting  and  processing  the 
information  to  reform  the  image,  nor  have  we  expanded  our  discussion  to 
include  the  orthogonal  (range)  dimension.  It  should  be  apparent,  however, 
that  the  synthetic  aperture  must  be  constructed  with  extreme  care  to  realize 
its  potential  fully.  For  example,  the  azimuth  resolution  suffers  if  the 
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satellite  is  perturbed  from  the  perfect  trajectory  by  a  significant  fraction 
of  its  operating  wavelength  (23  cm  for  Seasat)  as  it  traces  out  the  16-km 
aperture.  Various  types  of  aberrations  can  occur,  all  of  which  lead 
ultimately  to  loss  of  resolution  and  contrast  in  the  image. 

Having  explored  some  of  the  fundamental  requirements  for  aperture 
synthesis,  we  shall  now  summarize  systematically  the  major  steps  in  forming  an 
image  in  terms  of  a  stationary  point-source  response.  Because  an  imaging 
system  must  be  essentially  linear,  superposition  arguments  can  extend  the 
results  to  an  arbitrary  distribution  of  radar  backscatter.  A  point  source 
(step  1)  having  been  illuminated  by  a  coherent  radar,  emits  a  series  of 
concentric  wavefronts  (step  2).  Of  course,  the  emission  occurs  only  while  the 
point  source  is  within  the  aperture  beam,  as  discussed  above.  The  spacecraft 
cuts  through  the  concentric  wavefronts  (step  3)  and  the  SAR  receiver 
intercepts  an  energy  flux  that  varies  with  position  (or  time)  as  the 
wavefronts  are  traversed  (step  4).  This  wavefront  record  is  usually  referred 
to  as  the  Doppler — or  signal — history,  but  it  can  also  be  considered  a 
hologram,  diffraction  pattern,  or  a  one-dimensional  zone  plate.  The  essential 
point  is  that  this  wavefront  record,  which  contains  a  complete  phase  and 
amplitude  history  of  the  point  source  for  the  entire  synthetic  aperture 
interval,  also  contains  adequate  information  to  reproduce  a 
diffraction-limited  version  of  the  original  point  source.  The  wavefront 
record  is  transferred  via  a  data  link  (step  5)  from  the  spacecraft  to  any  of 
several  ground  stations,  where  it  is  recorded  on  either  digital  tape  or 
optical  film.  For  Seasat,  the  signal  was  normally  recorded  on  tape  at  the 
station  and  was  later  transferred  to  film  at  a  central  facility.  To  reproduce 
the  point  source  from  the  wavefront  record,  a  coherent  reference  function 
(such  as  a  family  of  plane  waves  from  a  laser)  (step  6)  impinges  on  the 


Doppler  history  (step  7).  The  Doppler  history,  which  is  essentially  the 
diffraction  pattern  of  a  point  source,  causes  a  lens-free  convergence  of  the 
plane  waves  in  free  space  (step  8)  and  forms  a  replica  of  the  original  point 
source . 

In  practice,  there  are  several  variations  to  the  simple  scheme  described 
above,  some  of  which  are  helpful,  some  nuisances,  and  others  inevitably 
destructive.  A  simple  lens  inserted  at  step  8  can  shorten  the  convergence 
distance.  On  the  other  hand,  the  spacecraft  trajectory  (step  3)  is  better 
described  by  an  arc  of  varying  center  of  curvature  than  by  a  straight  line. 
This  is  clearly  a  nuisance,  entailing  an  adaptive  processing  strategy. 
Atmospheric  turbulence,  reference  function  instabilities,  and  lens  aberrations 
are  examples  of  destructive  and  uncorrectable  sources  of  contamination.  Much 
of  the  effort  and  expense  of  spaceborne  SAR  can  be  attributed  to  the  need  to 
account  precisely  for  the  systemic  sources  of  contamination  and  to  minimize 
the  random  sources. 

The  extension  of  SAR  image  formation  to  the  range  (cross-velocity) 
direction  is  straightforward,  but  it  places  strict  timing  and  synchronization 
requirements  on  the  design  and  severely  restricts  the  total  range  interval  (or 
corresponding  ground  swath  width)  that  can  be  accommodated.  Range  information 
is  possible  in  a  SAR  only  because  the  synthetic  aperture  need  not  be 
continuous,  but  may  be  constructed  with  samples;  that  is,  the  transmitter  may 
be  pulsed.  It  is  enough  that  one  pulse  be  transmitted  each  time  the  real 
aperture  moves  by  half  its  length.  This  is  called  "filling  the  aperture"  and 
leads  to  a  maximum  time  interval  x  of  D/2v,  during  which  range  information  can 
be  collected.  For  D  =  12  m  and  v  =  8  km/s,  t  =  750  ps.  In  actual  practice, 
the  aperture  is  slightly  "overfilled"  to  reduce  the  possibility  of  spurious 
signals  that  could  lower  image  contrast.  The  Seasat  SAR,  for  example. 


typically  operated  at  an  interpulse  period  of  600  |Js.  Moreover,  somewhat  less 
than  half  of  this  time  interval  represented  signals  of  sufficient  quality  to 
produce  good  imagery. 

The  synchronization  constraints  imposed  by  the  Seasat  geometry  were  as 
follows.  Transmitter  pulses  are  emitted  every  600  |Js  in  the  cross~track 
direction,  20°  away  from  nadir.  The  pulses  form  concentric  expanding  rings  of 
energy  with  a  separation  of  about  200  km  (cT  =  3*10°  m/s  x  600*10  s  - 

200  km).  With  the  satellite  orbiting  at  an  altitude  of  800  km  there  are  four 

such  pulses  descending  at  any  time.  Since  the  antenna  illuminates  only  that 
region  around  20°  from  nadir,  no  significant  backscatter  occurs  except  at  a 
slant  range  of  about  850  km.  This  is  no  accident,  since  the  round-trip 

distance  of  1700  km  must  be  chosen  to  allow  the  return  from  a  particular  pulse 

to  occur  exactly  between  two  subsequently  transmitted  pulses.  The  geometry 
and  pulse  interval  for  the  Seasat  SAR  are  chosen  to  allow  the  return  from  a 
given  pulse  to  occur  exactly  midway  between  the  eighth  and  ninth  subsequent 
pulses.  The  middle  200  ps  of  return,  represents  a  swath  width  of  about 
100  km  for  the  Seasat  geometry.  Large  swaths  inherently  require  long  real 
apertures  to  adequately  accommodate  all  ranges  and  simultaneously  fill  the 
synthetic  aperture  properly.  This  can  be  a  serious  limitation  in  spaceborne 
SAR. 
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